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ABSTRACT 


Underwater acoustic communications networks serve to pass data between submerged 
nodes and a command and control center ashore. Range data recorded incident to regular 
transmissions between nodes afford the acoustic network an additional use as a 
navigational aide or tracking system for a vehicle operating in the network domain. 
Previous studies developed algorithms to solve for the position of the mobile node using 
node-to-node range data. An alternative localization algorithm is proposed. This study 
implements the algorithms for experimentation in simulation and an actual maritime 
environment. Reference nodes are deployed to establish a Seaweb network in Del Monte 
Lake on the Naval Postgraduate School campus. Experimental range data are recorded 
and used to plot the path of a surface vehicle towing a network sensor node. GPS fixes 
are simultaneously recorded for comparison of the calculated track to the true track. This 
comparison throughout the experiment provides a more definitive measure of the 
localization algorithms’ perfonnance than was possible in previous studies. 
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EXECUTIVE SUMMARY 


Underwater acoustic modems arranged in a fixed network can be used as a means of 
navigation for an underwater mobile node operating in the network domain. Ranges 
between stationary nodes and the mobile node can be obtained as a by-product of 
through-water acoustic communications. Various algorithms have been developed to 
calculate the position of an underwater vehicle using node-to-node ranges. Testing these 
algorithms with a surface vehicle allows direct comparison of the calculated tracks to that 
provided by the Global Positioning System (GPS). 

Seaweb is a Navy-developed acoustic network technology that uses underwater 
acoustic modems to communicate information to an onshore command and control center 
via a gateway node. The gateway node can be a moored buoy or a surface vehicle. Range 
data measured incident to regular transmissions between nodes afford the Seaweb 
network an additional capability as a navigational aide or tracking system. 
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Figure 1. Schematic of deployed underwater nodes and towed mobile surface node in 

Lake Del Monte. 


Ranges are calculated with the round-trip propagation time of transmitted signals 
and the sound speed in the environment. At regular time intervals during the vehicle’s 
run, the mobile node transmits a broadcast ping. Upon receipt at each fixed node, that 
node returns an echo including the node’s specific address. The mobile node records the 
ping transit time and node address for each return echo. 
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The electro-acoustic transducers are omnidirectional. Assuming a two- 
dimensional system, for each range measurement, there is a range circle about the fixed 
node that describes the vehicle’s possible position. The exact position can be found at the 
intersections of overlapping range circles. Many sources of systematic error exist 
including those due to the two-dimensionality assumption, signal processing delays, 
variations in sound speed, and uncertainty in the fixed node locations. Overlapping range 
circles will not intersect at one point and will yield ambiguous solutions, as seen in 
Figure 2. The localization algorithms developed in previous theses use different methods 
to choose the solution that is most correct. 


Intersecting Range Circles (or Sample Vehicle Position Intersecting Range Circles (or Sample Vehicle Position 



Figure 2. Overlapping range circles from five submerged nodes, and a close-up view, 
for a given position of the surface node. The average position is calculated with a 
difference linearization algorithm for synthetic data. 

The localization algorithm developed by LT Michael Reed, USN uses a difference 
linearization method, inspired by the GPS implementation. This method linearizes a set 
of three range circle equations to find the two u nkn own x and y coordinates. Solutions are 
found for all combinations of three range circles and averaged to obtain a final solution. 

Pairwise algorithms use sets of two range circles to calculate solutions. The 
correct solutions must be determined from each pair of ambiguous solutions. Solutions 

xviii 


















are found for all combinations of two range circles and averaged to obtain a final 
solution. The benefits of each method are discussed. 

Synthetic data are created to test the accuracy and robustness of the localization 
algorithms. Assumptions are made for the mathematical models, including straight-line 
propagation of sound and a two-dimensional representation, which closely approximate 
the shallow lake. Random error is introduced to each range measurement to test the error 
tolerance of each algorithm. Synthetic data are generated for random network 
configurations and random vehicle paths within the bounds of the lake dimensions. 
Synthetic results are also analyzed to find preferred network geometries. 

The difference linearization (DL) algorithm performs better than the pairwise 
algorithms studied in previous theses, but encounters problems when little range data are 
available per fix due to fewer available combinations of range equations. The method 
uses matrix inversion and also yields errors for certain network geometries. 


Actual and Calculated Vehicle Path for Simulated Data 



Figure 3. Graphical representation of the performance of the filtered pairwise 
localization algorithm for a sample synthetic track. 

Because of the shortcomings of the existing algorithms, a new “filtered pairwise” 
localization algorithm is proposed. The filtered pairwise (FP) algorithm uses all available 
combinations of two range circles to find a position solution. The method of calculating 
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the solution for two range circles depends on the particular geometry. In cases where the 
two range circles do not overlap, the point between the two range circles is chosen. In 
cases where two solutions exist, all solutions are compared to one another. The error 
between solutions is minimized to eliminate one of the two ambiguous solutions. A 
simple average of the remaining solutions is done to find the final estimated position. 

The algorithm is not mathematically complicated, but uses multiple logic 
structures for solution comparisons. The computation time is less than that for the 
difference linearization method. For simulations with ample range data per fix, the 
average error for the filtered pairwise method is larger. For situations with limited range 
data per fix, error due to matrix inversion in the difference linearization method makes 
the filtered pairwise method a good alternative. 

A Seaweb acoustic network was established in the Del Monte Lake on the Naval 
Postgraduate School campus in April 2012. The Seaweb network enabled testing with 
experimental controls. A series of experiments was performed with a towed surface node. 
The quiet boat towed a sled fitted with the mobile node and GPS loggers to 
simultaneously record GPS fixes and range data from the submerged Seaweb nodes. 
Localization performance with experimental range data can be quantified by comparison 
to the GPS track. 

Previous studies have shown Seaweb positioning to outperform inertial 
navigation. Submarines and unmanned underwater vehicles must surface periodically to 
reestablish their position via GPS due to errors inherent in inertial navigation. The goal of 
this research is increased accuracy of Seaweb range data and tracking algorithms to 
improve navigation of submerged vehicles. This capability can enable future deployable 
underwater ranges. 


xx 
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I. INTRODUCTION 


A. OVERVIEW 

Networks of underwater acoustic sensors are being advanced for a variety of 
applications including passive environmental data collection and surveillance. Networked 
underwater sensors can further be used as navigational aids for surface and subsurface 
vehicles. 

U.S. Navy Seaweb networks use fixed and mobile distributed nodes to transport 
communications to and from onshore coimnand and control centers. Seaweb modem 
firmware is used with commercial off-the-shelf acoustic modems that are adaptable for 
shallow or deep ocean operations, with use on surface vehicles, submarines and 
unmanned underwater vehicles (UUVs). The standard commercial modem firmware 
supports simple node-to-node communications. The Seaweb firmware supports network 
architectures and allows sensors to act as repeater nodes with link-layer and network 
layer protocols, described in Chapter II [1]. A description of Seaweb and references to 
past Seaweb operations are included in [2]. 

Seaweb modems can be anchored and suspended above the ocean floor, towed by 
a surface vehicle, fixed to a buoy, or integrated onto a mobile underwater vehicle. Each 
modem has a specific address that is included in all data packet transmissions. Once 
deployed, a network discovery is initialized that autonomously establishes an ad hoc 
network topology for communications. Anchored modems, referred to as repeater nodes, 
repeat infonnation through the network to a gateway node on a moored buoy or surface 
vehicle. The gateway node sends or receives information from a local or onshore 
operator. All network nodes can be accessed via the gateway node. 

In this thesis research, anchored modems communicate acoustically with a 
modem on a towed surface vehicle that is connected via serial port to a computer. 
Anchored modems are also referred to as remote nodes. The towed modem is referred to 
as the local node. 
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B. OBJECTIVE 

Seaweb communications provide node-to-node range information. A mobile 
vehicle equipped with a network node can utilize the range information from fixed nodes 
to track its position, using a range-based localization algorithm. Previous studies by [3]- 
[5] have shown this method to outperform the inertial navigation of an unmanned 
underwater vehicle and to avoid the need of periodically surfacing to reestablish its 
position via GPS, due to fix expansion inherent in inertial navigation. Increased accuracy 
of Seaweb range data and tracking algorithms will aid navigation of manned and 
unmanned submerged vehicles and can be the basis for future deployable underwater 
ranges. 

C. APPROACH 

This thesis reviews existing range-based localization algorithms. Previous thesis 
work by [3]—[5] included development and testing of localization algorithms with 
simulated and experimental node-to-node range data. The best-performing algorithm 
developed in [5] is tested with synthetic data and a new algorithm is proposed. 

A Seaweb network was established in the Del Monte Lake on the Naval 
Postgraduate School campus in April 2012. The Seaweb network enabled testing with 
experimental controls. A series of experiments were performed with a mobile surface 
node. A quiet boat towed an instrumented sled, fitted with the mobile node and GPS 
loggers. Ranges between underwater modems and the towed surface node were recorded 
while simultaneously recording GPS fixes. The experimental data are used to test the 
existing and new algorithms. The GPS track is a check on the performance of the 
localization algorithms. 
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II. ACOUSTIC RANGING 


Underwater ranging with acoustic signals is a challenge due to dynamic 
environmental conditions. Spatial and temporal variation of the medium influences sound 
propagation characteristics as described in [6], Further, anchored references are not 
stationary due to the dynamic environment. These effects and basic Seaweb acoustic 
communications are described. 

A. UNDERWATER SOUND PROPAGATION 

1. Sound Speed 

The sound speed in water is five orders of magnitude less than the speed of 
electromagnetic communications in air. This limits the bit rate of communications. Long 
round-trip propagation delays also introduce error for mobile sensors, whose movement 
during the course of ranging may be significant. Sound speed in water is highly variable 
and depends on salinity, temperature and depth. The surface layer is more dynamic than 
deeper ocean layers and will have greater temporal and spatial variations in sound speed. 

2. Geometric Spreading 

Sound energy spreads as it travels out from a source. This effect only depends on 
range from the source and is independent of frequency. When an acoustic wavefront 
spreads spherically from a source, the acoustic intensity is proportional to Mr , where r is 
the range from the source. Cylindrical spreading occurs at ranges where the sound is 
trapped in a channel by surface and bottom boundaries. The intensity amplitude of the 
wavefront is proportional to Mr. Geometric spreading limits the range of acoustic 
communications for a given source strength and receiver sensitivity. 

3. Attenuation 

In addition to the geometric drop in signal strength, signals are attenuated over a 
distance due to effects such as scattering and absorption. Attenuation depends on a 
number of variables including range, frequency, temperature, salinity, pH, viscosity, 
depth, and bottom and surface roughness. 
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4. 


Noise 


Man-made and ambient noise can limit the range of acoustic communications. 
Noise sources include machinery, cavitation, shipping, flow noise, electronic noise, noise 
due to the sea state, and biologies. Noise levels vary for different frequency bands. 

5. Refraction 

Waves that are obliquely incident on a boundary between two media will bend 
according to Snell’s law 

c, c 9 

— 1 — = —-— = constant (l) 

cos 6 ] cos 0 2 

where c is the sound speed and 6 is the grazing angle of the ray measured from 
horizontal. For a constant sound speed gradient, waves will follow circular paths with a 
radius of 

R _ c O) 

\g\ cos 0{z) 

where z is the depth, c(z) is the depth-dependent sound speed, g is the sound speed 
gradient defined as g = Ac/Az and 0 is the angle from horizontal at the depth z. Typical 
radii of curvature for sound rays are on the order of kilometers. Thus, for short ranges, 
straight line propagation of acoustic signals is a good estimate. 

6. Multipath 

Multipath signals occur when the transmitted signal follows multiple paths to the 
receiver, each with a different propagation distance and travel time. When using 
propagation time to determine range, multipath can lead to over- or underestimations. 
Multipath effects are described in [6]-[8], The phenomenon occurs due to the omni¬ 
directionality of the transducer and boundary reflections. If the source or receiver is near 
the bottom or surface boundary, reflections arrive near in time to the direct path signal 
and can cause destructive interference. Impulse responses show that the amplitude of 
reflected signals can sometimes be stronger than the direct path signal. Multipath signals 
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can also arrive earlier than direct path signals. If obstacles obstruct the line of sight from 
source to receiver, multipath signals can be mistaken for the direct signal. 

Increasing the distance of modems from the bottom and surface boundaries 
usually will decrease multipath interference. Intersymbol interference (ISI) can be 
minimized by lowering the data rate and increasing symbol duration. Guard times 
between symbols allow multipath signals to die out. 


Doppler Spreading 


Doppler spreading occurs when the source and receiver are in motion relative to 
one another. A given frequency is shifted from the source frequency,/,, according to 


f = fo 


1 +’ 


1 -’ 


( 3 ) 


c J 


where v r is the receiver velocity, v s is the source velocity, and c is the sound speed in the 
medium. Multipath signals that arrive at the receiver at different angles will experience 
different Doppler shifts due to the projection of the relative motion. Different frequencies 
also experience a different Doppler shift. 


B. FURTHER ENVIRONMENTAL CONSIDERATIONS 
1. Location of Deployed Nodes 

Difficulties in node deployment cause uncertainty about the node’s position on 
the sea floor. [9] discusses the trajectory of a deployed node as it travels to the sea floor 
when ocean currents are present. 
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Figure 1. Pictured is a model trajectory of a deployed sensor. The graph shows the 
average displacement of sensors from computer simulations as the velocity of the 
current is increased. The trajectories of two types of sensors, labeled “Sensor” and 
“Uw gateway,” are modeled at three bottom depths. From [9]. 


The trajectory equations assume complete knowledge of the ocean currents as 
they vary with depth. This knowledge of deployed nodes is impractical, but simulations 
give a good estimation of the horizontal displacement to be expected. Shallow sensor 
deployments are much less impacted by currents. The difference in “Sensor” and “Uw 
gateway” trajectories, shown in Figure 1, vary based on the shape and weight of each 
particular object. 

2. Node Mobility 

When anchored, nodes suspended above the ocean floor will drift due to currents. 
Buoys can also be affected by wind forces. Tides and currents can affect the location and 
orientation of modems and surface floats as shown in Figure 2. The area in which the 
modem might be located is bounded by a watch circle. The radius of the circle is 
determined by length of the tether and the water depth. The anchor position may also 
migrate due to a changing sea floor. A modem on a surface buoy can use GPS to update 
its position, but the exact position of submerged nodes is difficult to detennine. 
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Figure 2. Tides and currents will affect the location and orientation of an anchored 
modem. The modem can lie anywhere within a watch circle, whose radius 
depends on the length of tether and the water depth. 

C. SEAWEB RANGE DATA 

1. Format 

The local node can send a number of 9-byte specialized data packets, termed 
utility packets. In this thesis, the basic ping command is used, which first sends a ranging 
signal from the local node to a remote node, followed by a PING utility packet. A 
matched filter at the remote node detects the ranging signal, which is a Flyperbolic 
Frequency Modulated (FIFM) chirp. The time of arrival at the remote node is identified at 
the peak amplitude of the incoming chirp. In most cases, the multipath signal amplitudes 
will be lower than the direct path signal and will not lead to false ranges. 

The remote node replies with a ranging signal and ECHO utility packet. The local 
node uses the time difference between the outgoing signal and the incoming signal to 
calculate the round-trip propagation time, termed the “delay” time, and range. The delay 
time already accounts for the processing time at the remote node. Procomm, by 
Symantec, is a terminal emulation tool that can be used for data communications, as 
described in [10]. Procomm is used to send commands to the local modem and display 
the received information. A sample ping command and response displayed in the 
Procomm terminal are shown in Figure 3. 
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atep20 

Send Ping-time out: 009.9 secs user:ll> 

XMT:[20]<[26] PING Time: 19:57:45.6796 
RCV:[26]<[20] ECHO Time: 19:57:48.2218 
LOCAL CMD RESPONSE (63) : 

Range Data: Node 26 to Node 20; Delay 151.0 ms; Range 113.2 m 


Figure 3. Procomm terminal display of information from Seaweb utility packet. The 
atep command sends a ping to Node 20 from Node 26. The “Delay” is the round- 
trip propagation time for the transmission. 


The range value is calculated with the delay time and a sound speed of 1500 m/s. 
This is intended for use in an ocean environment, but provides an upper estimate of the 
range when operating in freshwater. The delay time value is given to a resolution of 0.1 
milliseconds. Immediate subsequent range measurements, taken between two submerged 
modems that are approximately stationary, show variation of at most 1 millisecond in the 
delay times. For a sound speed of 1481 m/s, 1 millisecond error corresponds to a range of 
1.481 meters. 

2. Time Synchronization 

For certain localization schemes, time synchronization among the network nodes 
is critical. The accuracy of clocks will drift over time. The clock times can vary by an 
offset or can be skewed by a scalar factor. 

When the round-trip delay time is calculated by using both the arrival time at the 
remote node and then the arrival time at the local node, error from clock differences can 
be introduced. However, if the round-trip delay time is calculated at the local node using 
transmission and arrival times, as is done in all methods considered here, no clock 
synchronization is needed. 

3. Basic Range Calculation 

The range from the remote node to the local node can be estimated as 

r ij = C ^ t ab/ 2 ^ ^ 

where r tj is the range between the local node i and the remote node j, c is the sound speed 
for the environment, and A t a b is the round-trip delay time calculated at the local node. 
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The Broadcast Ping Seaweb command transmits a signal from the local node to 
neighboring remote nodes in the network. The local node then waits for responses for a 
set period of time. The remote nodes each have an assigned random dwell time. After 
receiving a broadcast ping, the remote node waits for its specific dwell time before 
replying to the local node. The dwell time is used to ensure return signals from the 
neighboring nodes do not interfere with one another. The dwell time is subtracted from 
the elapsed time at the local node, which then computes the range. 


nodei nodej 



RTS/CTS Handshaking 


nodei nodej 



Utility Packet sent 
with Broadcast Ping 

Handshaking disabled 


Figure 4. Seaweb handshake operation and broadcast ping with dwell time. After [3] 

and [5]. 


Using Equation (4) to calculate ranges requires a number of assumptions. The 
equation is true if the speed of sound is spatially and temporally constant, sound 
propagates in a straight line, the remote nodes and the local node are stationary, and there 
is no error in the calculation of the elapsed time. 

The Seaweb firmware uses handshaking when transmitting data. The local node 
will send a request-to-send (RTS) message to a remote node that takes the remote node 
out of its idle, low-power state. The remote node replies with a clear-to-send (CTS) 
message to indicate it is ready to process incoming data. The local node then sends the 
data packet. If errors exist in the data packet, the remote node sends a selective automatic 
repeat request (SRQ) for the portion of data with errors. The local node replies with the 
truncated data. This SRQ/DATA exchange repeats for a set number of times. 


9 














When handshaking is used in operational situations, range data are collected for 
all RTS/CTS transmissions. For the experiment discussed in this thesis, all range data 
were collected using the PING/ECHO dialogs rather than from RTS/CTS handshaking. 
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III. LOCALIZATION TECHNIQUES 


A. RANGE-BASED APPROACHES 

A variety of localization schemes exist for through-air communications, which do 
not directly translate to the underwater environment. Sensor networks and localization 
schemes must be designed to work in a spatially and temporally changing medium. Most 
underwater acoustic localization techniques use range measurements with static 
references. A 2010 survey of localization schemes by [8] and [11] classify range-based 
approaches into three categories. 

1. Received Signal Strength Indicator (RSSI) 

The RSSI approach to localization uses the strength of incoming signals to 
determine the distance to a source. The source strength must be known and a reliable 
model for the attenuation over range is needed. The RSSI method for range- 
determination is not ideal in an underwater environment due to the many variable 
contributions to transmission loss and the existence of multipath 

2. Time Difference of Arrival (TDOA) 

The TDOA approach is commonly used in localization problems. This approach 
requires two separate time-delay measurements between two objects to calculate a range. 
For example, a node might send an radio transmission and an acoustic signal 
simultaneously. The range to a receiver is calculated based on the known propagation 
speeds of the two signals and the difference in arrival times. A related technique is that 
described in [12]. Distances to an object are found with the time delay between incoming 
multipath signals from the same transmission. This method requires a good model of the 
multipath propagation and the geometry of the nodes. In general, TDOA approaches are 
not well suited for a variable underwater environment. 


11 



3. Time of Arrival (TOA) 

The commonly used scheme for underwater localization is the TOA approach. 
This uses transmission times between nodes and the speed of propagation to calculate a 
distance. If the ranges are calculated using one-way transmissions, clock synchronization 
is required. The accuracy of this approach requires accurate knowledge of the sound 
speed and the length of sound propagation paths. The data needed for the TOA approach 
is already collected with Seaweb sensor network communications, thus the TOA 
approach is used for thesis. 

B. TOA APPROACHES TO LOCALIZATION 

1. Intersecting Circles and Spheres 

When using omnidirectional transducers, the direction of the remote node with 
respect to the local node is unknown. The range calculated with Equation (4) is the radius 
of a sphere centered at the remote node, on which the local node is located. In two 
dimensions, the range is the radius of a circle, herein called a range circle, centered at the 
remote node. For localization with Seaweb range data, the positions of the remote nodes 
are assumed to be known. Overlapping range circles will pin down the location of the 
local node. In general, three overlapping circles are needed to uniquely find the local 
node. For three dimensions, four overlapping range spheres are needed. The solutions for 
overlapping range spheres are worked out in [13]. Solutions for overlapping range circles 
are found in Chapter IV. 

2. Intersecting Hyperbolas and Hyperboloids 

An alternative method to using overlapping circles is described in [11]. If the 
range measured by a vehicle to the node N\ is R] and the range measured to node Ni is Ri_, 
the vehicle will be located somewhere on a hyperbola for which the value {R\-Rt) is 
constant. [11] proposes using intersecting hyperbolas, or hyperboloids in three 
dimensions, to find the vehicle’s position. Two intersecting hyperbolas will always yield 
one solution, versus the intersecting circles approach that gives two ambiguous solutions. 
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3. Unambiguous Solutions 

Three methods for finding an unambiguous solution for overlapping range circles 
are discussed in Chapter V. A fourth, commonly used approach is the Least Squares (LS) 
algorithm, which gives best position estimate by finding the solution that minimizes the 
squared range error. The LS approach is used to solve non-linear equations and to 
circumvent the problem of ambiguous or no solutions when using overlapping range 
circles. An LS algorithm is also used for minimizing error with intersecting hyperbolas. 
The range error, s h for node i can be defined as 

A = n ~ yjix-xf +(y-y i ) 2 (5) 

where (x, y) is the solution, (x h y t ) is the remote node location, and r, is the measured 
range. [14] describes various methods of solving for 

N 

min^r. (6) 

/=1 

Typically, the range equations are approximately linearized and solved iteratively or 
directly with matrix calculations. As described in [15], the LS algorithms are good for 
cases when range error values cannot be assigned a confidence value or weight based on 
a predictable error model. [15] investigates the sensitivity of the intersecting hyperboloids 
and intersecting spheres localization approaches with LS algorithms. [11] proposes an 
alternative to the LS method for intersecting hyperbolas, however [15] and [13] find the 
intersecting spheres method is much more tolerant to errors in the range data in general. 

For this thesis, the localization problem is two dimensional. Intersecting circles 
are used. An LS algorithm is not used. The handling of ambiguous solutions and non¬ 
intersecting circles for sets of two range circles is described in Chapter IV and 
Chapter VI. The difference linearization method is also investigated and described in 
Chapter V. 
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IV. VEHICLE TRILATERATION WITH RANGE DATA 


Ideal range data will yield range circles that overlap at a point. In practice, errors 
in range data will yield circles that overlap at two positions or range circles that do not 
overlap. The solutions for each case are calculated and related errors are discussed. 

A. GEOMETRY 

The intersection points of two overlapping circles with known center positions 
and known radii can be calculated as shown. 



Figure 5. Overlapping range circles can yield two ambiguous solutions, labeled here as 
P2. The positions of Node i and Node j are known as well as the radii. After [16]. 


Point P 1 is first defined with coordinates {xp\, y P i) between the two nodes as 
shown in Figure 5. The distance from PI to Node i, with coordinates (x h v ; ), is defined as 
a. Likewise the distance from P 1 to Node j, with coordinates (xj, yj), is defined as b. The 
intersection points are a distance h from point P\ as shown. This leads to the 
relationships 

R 2 = a 2 + h 2 

R 2 =b 2 +h 2 

J 

and 

d = a + b = V(Ax,) 2 + (A y.j) 2 (7) 

where Axy and zhy are defined by 

A x ■ ■ = x ■ — x ■ 

IJ l J 
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V'\ ; = .'■/ .'■/• 

First solve for the unknown variables a and h. 


(d-af =b 2 


d 2 -2ad + a 2 =R 2 -h 2 
d 2 -2 ad + 0 1 = R 2 -R i 2 + 


a = 


(R 2 -R 2 +d 2 ) 


2d 


><=d, 


2 - a 2 


The center point can now be found using the relationship 

(/j - Node t ) _ a 

( Node. - Node j ) d 


from which it can be shown 


a 


x Pi= x i + j( x j- x i) 


y P i = yi+j(yj-yi)- 


( 8 ) 


A relationship for the perpendicular lines in terms of P2 and known variables is written 

(Xp 2 ~ Xp\ ) _ — (Vp 2 — Tpj) _ i h 


to find 


O ,-yd ( x r x i) 


Xp2=x Pl ±^(y j -y i ) 

yp2=y n +^( x j- x i) 


d 


( 9 ) 

( 10 ) 


This fonnulation from [16] is general for any choice of Node i and Node j. 


If the circles do not intersect, the value of h, calculated with Equation (8) will be 
imaginary. For such cases, a variable c is defined to be the distance between the two 
circles. Using the distance d calculated in Equation (7), c is 

c = d - R . - R j 

If c is greater or equal to zero, the circles are oriented as shown in Figure 6. 
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Figure 6. Range circles may not intersect. In the case shown, the point between the two 

range circles is chosen as the solution. 


In this case, the solution is the point that is located halfway between the circles. 
The distance a is now defined as 


a = II +(c/ 2). 


The relationship 

a _ (*pi -Xj) _ ( y n -yj) 

d (Xj -x) (y. -y.) 

is used to find 

= x i + -(x j -x i ) 
d 

ypi=yi + ^(yj-yi) 

d 


( 11 ) 


( 12 ) 

(13) 



Figure 7. Range circles can lie within one another as shown. The point at which the 

circles are closest is chosen as the solution. 
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Cases where h is imaginary and c is less than zero correspond to two non¬ 
intersecting circles that are inside one another. The solution is located at the point 
halfway between the circles where the circles are closest. Assuming R ; is greater than R h 
the distance e is defined as 

e = R t — R . — d 

from which it is found 

a = R j + (e/2) 

Xpi = *i + “7 (*/“*/) (14) 

d 

Vpi = yj + jiyj-Vi) • (15) 

B. INHERENT AND SYSTEMATIC ERRORS 
1. Geometric Dilution of Precision 

Small errors in just one range can quickly degrade the accuracy of these 
calculations as shown in Figure 8. Geometric dilution of precision (GDOP) refers to the 
geometry-sensitive change in error. The position solutions from certain sensor network 
geometries are less sensitive to error than others. Thus, it is difficult to calculate a 
confidence value for each solution. 




Figure 8. An example of GDOP. A small error in one range measurement will lead to a 
large error in the calculated solutions for certain geometries. 


2. Two-dimensional Model 

All algorithms used in this thesis assume a two-dimensional system. This 
approximation is valid in the shallow lake. Fathometer measurements showed the lake to 
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be between 4 and 8 feet deep. The mobile local node was suspended 13 inches beneath 
the lake surface. The deepest node was suspended 48 inches below the surface. Thus, if 
the boat is maneuvered over the deepest node, the range measurements have an error of 
35 inches. All other errors due to the two-dimensional estimation are less than this. 

Because straight-line propagation of sound is assumed, this error is corrected by 
projecting all ranges onto a two-dimensional plane. The converted ranges are used as 
inputs to the localization algorithms. 



Figure 9. Projection of three-dimensional range onto two-dimensional coordinate 

system. 


A range between nodes i and j, defined as D , has magnitude 

r iJ3D =^Ax ij 2 + Ay i / + Az i / (16) 

The projection of the range onto the x-y plane is defined as d and is shown in Figure 9. 
We can rearrange Equation (16) to solve for the projected range 

2 2,2 2 2 , , 
V ij2D ~ X ij + ) ; ij ~ r iflD ~ Z ij (17) 

r ij2D=J r j3D- Z i/ _ ( 18 ) 

This correction is used, but makes a minimal difference in calculating the path of the 

vehicle for such a shallow geometry. 

3. Vehicle Motion 

If the broadcast ping is used to communicate with N remote nodes, the remote 
nodes have no built-in dwell times, and the vehicle is stationary, the echoes will return to 

the local node at different times due to the differing ranges. If the vehicle is moving, the 
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vehicle will receive the N echoes at N different positions along its track. The return 
propagation times will change, and will be reflected in the round-trip propagation time. 
The time delay data received at the N different positions are used to calculate one 
waypoint that corresponds to the time the broadcast ping was sent. Thus, the vehicle’s 
motion introduces two separate errors to the calculation of position: the ranges are 
inexact, and they do not correspond to the same position and time. 


outgoing ping travels At x 

incoming echo travels At 2 


<=> X □ 

Xi x 2 

-► 

vehicle motion from x 2 to x 2 
over t = At 2 + At 2 — 


X 


-A 

fixed node at x N 


Position calculated 
due to vehicle motion 


Figure 10. Schematic of the effect of vehicle motion on ranging. 


The most simple and dramatic example of the first error occurs when the vehicle 
is moving directly toward or away from a fixed node. As pictured in Figure 10, the range 
calculated corresponds to neither the position when the ping was sent, nor the position 
when the echo was received. If the echo is sent the moment the ping arrives at the remote 
node, these errors are usually negligible because the propagation speed of the signal is 
much greater than the vehicle speed. Note that Figure 10 is not to scale. 


outgoing ping travels Atj 

incoming echo travels At 2 

incoming echo travels fits 

O Q D X_CD_A fixed node at Xn 

xi x 2 xs with dwell time, t 

-► vehicle motion over t = Ati + At 2 , x = 0 

-*■ vehicle motion over t = Ati + x + At,, x * 0 


o 


Position calculated 
due to vehicle motion 


x 


Position calculated due to 
vehicle motion and dwell time 


Figure 11. Schematic showing the effect of built-in dwell times for the remote nodes. 

The addition of a dwell time yields a range error equivalent to that found when 

the vehicle speed is increased. 
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When a Broadcast Ping command is used, the remote nodes send a reply after a 
built-in dwell time, r. This dwell time increases the ranging error as shown in Figure 11. 
Processing time at the local and remote nodes increases the ranging error in the same 
way, as does increasing the vehicle speed. All of these variables should be minimized for 
experimentation. 
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V. PREVIOUS WORK ON VEHICLE LOCALIZATION 

ALGORITHMS 


A. WEIGHTING METHOD 

The algorithm developed by [3] uses a weighting method to estimate the position 
of the local node. Each possible solution from two overlapping range circles is assigned a 
weight based on the proximity of all other solutions. Outlier solutions are given smaller 
weights. Clustered solutions are assigned higher weights. The solutions are multiplied by 
the weights, summed, and divided by the sum of the weights to find the estimated 
position. 



In Equation (19), (x, y ) is the estimated position, the (xj, V/)’s represent all the raw 
solutions, and the W/s are the calculated weights. 

B. CENTER OF MASS METHOD 

The center-of-mass method developed by [4] calculates the vehicle position by 
averaging all the solution values to find the center of mass for the system. Solutions 
distant from the center of mass are discarded and the center of mass is recalculated. 

C. DIFFERENCE LINEARIZATION METHOD 

The difference linearization (DL) method, adapted from a GPS localization 
method, is investigated in [5] and found to perform better than the weighting or center of 
mass algorithms. The algorithm uses combinations of three range circles to find a 
solution. Three range equations can be written as 

r 22 f(x-x t ) 2 +(y~y 1 ) 2 ^ (x 2 - 2xXj + x t 2 + y 2 -2yy x + y 2 
r{ = (x-x 2 ) 2 +(y~y 2 ) 2 = x 2 -2xx 2 +x 2 2 +y 2 -2yy 2 +y 2 2 (20) 

r 2 ) [(x-x 3 ) 2 +(y-.y 3 ) 2 J [x 2 -2xx 3 +x 2 +y 2 -2yy 3 +y 2 ^ 

Subtracting the second range equation from the first gives 

r { 2 - r 2 2 = - 2 x x t + x { 2 - 2 y y x + y { 2 + 2 x x 2 - x 2 2 + 2. y y 2 - y , 2 
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which is linear in x and y. A pair of simultaneous linear equations gives the matrix 
equation 


( 2 

2 > 

r i 

~ r 2 

2 

2 

VC 

-C J 


= 2 


(* 2 “*i) (y 2 -y l ) 


\ / \ f 2 2 , 2 

1 x l -x 2 +y x -y 2 


v (x 3 -x 2 ) (y 3 -y 2 ) 

that has an exact solution of 


x 

w 


+ 


2 A 


z z . z z 

x 2 -x 3 + y 2 -y 3 


fX) 

_ 1 

'(x 2 -x l ) (y 2 -y l )^ 

-1 

( 2 2,2 2,2 2 \ 
r \ ~ r 2 + x 2 — Xj +y 2 y x 


~2 

V (x 3 -X 2 ) (y 3 -y 2 )j 


2 2 2 2 2 2 
\ r 2 x 3 — x 2 +y 3 ~y 2 j 


( 21 ) 


( 22 ) 


This can be written as 


R=-A^C 

2 


(23) 


where R 


v 


It; 


A= 


(x 2 -^q) (v,-v,)' 
^(x 3 -x 2 ) (y 3 -y 2 ) 


, and C = 


( 2 2,2 2,2 2 \ 

C + * 2 +T 2 -Ti 

2 2,2 2,2 2 

r i ~r 3 + x 3 -x 2 + y 3 -y 2 


This method also works for overdetermined cases where there are more than three 
ranges available. If there are N fixed nodes, the matrix A will have size (N -1 x 2) and the 
vector C will have (N- 1) elements. It was found in [5] that calculating one solution with 
this method for more than three nodes was less accurate that averaging the solutions from 
all combinations of three nodes. Combinations of three nodes leave the matrix A a square 
matrix that can be inverted for an exact solution. 


1. Eigenvalues 

The DL algorithm runs into errors given certain configurations of the anchored 
modems. As seen in Equation (23), to calculate an x and y position, the matrix 


A= 


(x 2 -x l ) (v 2 -v,) 


must be inverted, and depends solely on the positions of three 


v (x 3 -x 2 ) (y 3 -y 2 l 

remote nodes in relation to one another. If the nodes are co-linear, the matrix will be 
singular with a zero eigenvalue, and cannot be inverted. If the nodes are arranged in 
approximately a straight line, the matrix will have a small eigenvalue. When solving 
linear equations with matrix coefficients, for data with non-zero error, a small eigenvalue 
will inflate that error. [17] shows that the inverse of a matrix can be represented by 
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(24) 



where M is an nxn matrix and e, is the eigenvector of the matrix corresponding to the 
eigenvalue A/, and T indicates the transpose of the vector. If the vector C from Equation 
(23) has error, this error is magnified with a small eigenvalue. In overdetermined 
situations, solutions can be filtered out that correspond to eigenvalues much less than one 
can be discarded before averaging. 


2. Condition Number of Matrix 

Errors also occur when the magnitude of the eigenvalues are well behaved. These 
errors can be avoided by checking the condition of matrix A. The condition number of a 
matrix is a measure of the accuracy of linear calculations using that matrix. The condition 
is a characteristic of the matrix itself, not of the calculation. Consider the calculation of 
solution vector x with the square coefficient matrix A and the constant vector b. 

A ~x = b (25) 

As described in [18], if there is a small error in the values of the coefficient matrix 
or in the constant vector, a well-conditioned matrix will have small errors in the solution 
vector. An ill-conditioned matrix yields large errors in x for small perturbations of the 
values in A or b. The condition number of a matrix is the same for the matrix inverse. 
The condition of an nxn matrix is defined as 

cond(A) = | A|| • || A -1 || > 1 (26) 

where | |A| | is the 2-norm, or Euclidean vector norm, of the matrix. As with vector norms, 
the norm of a matrix is a measure of the size of the matrix. The 2-norm is defined by [19] 
as 

Ax 

II All =max - 2 (27) 

2 x^O £ 

2 

or more simply 



9 . . j. j- 

where erf is the maximum eigenvalue of A A, A being the conjugate transpose of A. 
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Reference [18] proves via the definition and properties of norms that the error in 
vector x can be related to a perturbation in matrix A by 


Ax 


< 


x + Ax 
II ii 

and similarly for a perturbation of vector b 


(IN'II A- 


lllJ AA ll 


(29) 



Ax 


M 

J 

<1 

x + Ax 

-Mil ll 1 

\) 

b 


(30) 

This gives a mathematical representation of the condition number ^||a||-||a _ i ||), which is 
the upper limit of how error in one measurement scales error in x. 


[20] describes a rule of thumb to estimate error in a measurement. If m is defined 
as 

m = logio[cond(A)] (31) 

m represents the number of significant figures lost in the solution vector x due to a 
perturbation. Condition numbers close to unity are best. For overdetennined situations, 
the DL algorithm is edited to keep only those matrices with condition numbers less than 
10, corresponding to a one digit loss in accuracy. Matrices with low eigenvalues will also 
be ill-conditioned. Thus, the condition of a matrix can be used as a metric to identify 
good node geometries. The DL algorithm edited to check the A matrix condition number 
is referred to as the difference linearization with condition check (DLC) method. 


3. Number of Available Ranges 

A further limitation of the DL method is the need for three range circles to 
calculate one position. If the A matrix is ill-conditioned, it will be bad for all points on 
the track. Whereas in an over-determined situation, solutions from ill-conditioned 
matrices can be ignored, if only three remote nodes are available, all solutions will have 
inflated error. 
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VI. THE FILTERED PAIRWISE ALGORITHM 


A. ALGORITHM DESCRIPTION 


As an alternative to the difference linearization (DL) method and the previously 
developed pairwise algorithms, a fdtered pairwise (FP) algorithm was developed. With 
only three nodes available for ranging, pairwise algorithms generate three solutions from 
the three combinations of two range circles. Error is reduced in the pairwise algorithms 
by averaging the three solutions. The number of combinations available can be described 
with 


\ k j 


n\ 


(32) 


k\(n—k)\ 

read as ‘n choose k.' With five nodes, each method has an equal number of combinations 
of nodes. With six nodes, the difference linearization method has more combinations for 
averaging. In most underwater acoustic networks, six nodes will not be available for 
ranging. 


The pairwise methods do not encounter the complications of matrix inversion. 
The FP method uses simple geometry to calculate all solutions and uses logic loops to 
discard the erroneous solutions. 


B. ALGORITHM STRUCTURE 

For each waypoint, an estimated position is calculated from all combinations of 
range circles as described. 


1. Computing All Solutions 

Solutions for range circles that do not intersect are calculated first, as shown in 
Chapter IV. These combinations have only one solution that, even with GDOP error, is a 
good initial guess at the vehicle position. The solutions are saved. Next, for all cases with 
intersecting range circles, the two ambiguous solutions are calculated and saved. 
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2. Comparing All Solutions 

It is assumed that for every pair of ambiguous solutions, one solution is ‘correct’ 
and the other solution is incorrect, and the correct solutions are near one another. For the 
first pair of ambiguous solutions, Si a and Sib, the distance between each solution and 
every other solution, correct or incorrect, is calculated. For example, given two 
ambiguous solutions for a different combination, S ia and Sw, the distances 

4=(^-SJ Q=(^-SJ D,=(S lh -S ih ) 

are calculated. One of the four distances represents the distance between two correct 
solutions. The smallest distance of A t and Bj is saved for point Si a , and the smallest 
distance of C, and D, is saved for Su,. This is repeated for all other pairs of ambiguous 
solutions and combinations of pairs. Finally, the saved distance values for S\ a are 
summed and compared to the sum of the distances for Sib- The smallest sum corresponds 
to the correct solution. This one correct solution is then used to choose the correct 
solutions from all other pairs of ambiguous solutions, using the distances already 
calculated. 


3. Averaging All Filtered Solutions 

The solutions from non-intersecting range circles and the chosen correct solutions 
are averaged to find a single position for the waypoint. No weights are applied, and at this 
point, further outliers are not discarded. 
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VII. ALGORITHM IMPLEMENTATION FOR SYNTHETIC DATA 

A. SETUP AND ASSUMPTIONS 

To test the algorithms, random tracks are created that are continuous and stay 
within given bounds. The dimensions of the simulation area are the approximate 
dimensions of the anticipated experimental area in a shallow lake. A set number of nodes 
are placed in the domain and are given random fixed positions. A mobile node issues 
broadcast pings at a given period along the track. The nodes reply with an echo after a set 
dwell time. The echo arrival times are used to calculate ranges. Random error or bias can 
be added to the range data. The range data are saved to be used as input to the algorithms. 
All variables can be edited as needed in the Matlab script. 

The simulations assume a completely two-dimensional scenario. There is some 
inherent error due to interpolation and limited resolution in Matlab when making the 
tracks and simulating the traveling pings and echoes. The errors are small and have a 
minimal effect on the perfonnance of the algorithms. 

B. CASE STUDIES 

Case studies comparing the filtered pairwise (FP) method, difference linearization 
(DL) method, and the difference linearization with condition check (DLC) method are 
run using five remote nodes. Again, with five nodes, there is an equal number of 
combinations of nodes, and thus, an equal number of possible solutions for each type of 
algorithm. A script is created that generates a random path with random node positions. 
Range data are generated and the average localization error for the algorithms is found by 
averaging the difference between the actual vehicle positions and calculated positions for 
the entire track. 100 iterations are performed for the script. 

1. Data Sets with Zero Range Error 

Simulations are first run with zero added range error and no dwell times for the 
nodes. Each algorithm is run with the same data. There is still error in the data due to the 
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motion of the vehicle. This error is minimal however, as the echoes are allowed to arrive 
at the same time and the pings take only milliseconds to and from the remote nodes. 


Localization Error for 100 Iterations of Error-Free Synthetic Data 


FP Method 

DL Method 

DLC Method 

Minimum Error (m) 

0.1305 

0.1001 

0.0926 

Maximum Error (m) 

0.7141 

4.8981 

0.3004 

Average Error (m) 

0.1978 

0.4743 

0.1458 


Table 1. Comparison of localization error for the three algorithms using error-free 

synthetic data. 


For error-free synthetic data, the FP method performs better than the DL method. 
The DLC algorithm has the lowest average localization error. The difference between the 
DL and DLC methods shows the effectiveness of checking the condition number of the 
inverted matrices. 

2. Data Sets with Random Range Error 

Error can scale with range or be an additive factor, depending on the source of the 
error. Ranging error can be due to indirect sound propagation paths, refracted 
propagation, motion of the remote nodes, motion of the vehicle, processing time at the 
remote or local node, and a number of other factors. Errors also occur due to interpolation 
and limited resolution in Matlab when creating synthetic data. Because the error is not 
predictable, random error is simply added to the calculated ranges. The rand command in 
Matlab is used to generate a random number between 0 and 1. The value is shifted and 
scaled to give a random number between -1 and 1. This number is then multiplied by a 
variable scalar factor. The scalar factor represents the maximum range error in meters. 
The equation used in Matlab and the equivalent expression are 

R e =^ o + 2£'(rand-0.5) 

\r e -r 0 \<e 

where Re is the range with error, R 0 is the original range, and E is the maximum added 
range error. Simulations are run for scalar factors of 0.5, 1, 1.5 and 2 meters. 
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Localization Error for 100 Iterations of Synthetic 


Data with Introduced Error 

45 -— - 



± 0.5 m ± 1.0 m ± 1.5 m ± 2.0 m 


Introduced Errorto Individual Range Measurements 


Figure 12. Comparison of localization error for the three algorithms using synthetic data 

with introduced error. 


When error is added to the synthetic range data, the FP and DLC methods have 
similar average localization errors that scale almost linearly with the maximum error 
added to individual range values. The improvement of the DLC over the DL method is 
evident. The FP appears to perform better with greater error. To investigate this, 10 
iterations of the algorithm are run over a range of introduced error values. 


Average Localization Error with Introduced 
Error for 10 Iterations of Synthetic Data 



Maximum Introduced Error to Range Measurements (m) 


Figure 13. A measure of the robustness of each algorithm is found by plotting the 
average localization error as it increases with the amount of error introduced to 
each individual range measurement. 
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The graph in Figure 13 gives a rough measure of the robustness of each algorithm. The 
average localization error for the FP method grows approximately linearly with the 
amount of error introduced to individual measurements. The growth of average error for 
the FP method is less than that for the DLC method, indicating that it may be the better 
choice for practical implementation. 

3. Data Sets with Dwell Time 

Data are created for which the remote nodes are assigned a random dwell time. 
The maximum dwell time for a node is set to be one half the time between pings, or 2.5 
seconds for these simulations. This is an arbitrary choice for the synthetic data. 


Localization Error for 100 Iterations of Synthetic Data with Dwell Times 


FP Method 

DL Method 

DLC Method 

Minimum Error (m) 

0.2596 

0.1956 

0.1744 

Maximum Error (m) 

1.0612 

55.6238 

0.8228 

Average Error (m) 

0.4802 

1.8873 

0.4542 


Table 2. Average localization error for the three methods when a dwell time is used for the 

remote nodes. 


The data have no artificially introduced error. The DLC method again performs 
best in this case. The effect of adding dwell times can be seen by comparing the error 
values with those from error-free range data in Table 1. 


Average Localization Error with Dwell Times 
using 10 Iterations of Error-free Synthetic Data 

1.2 



0 1 2 3 4 s 6 

Maximum Dwell Time for Remote Nodes (s) 

Figure 14. The performance of each algorithm with dwell times at the remote nodes is 

plotted. The synthetic data are created for a vehicle moving 0.5 m/s. 

32 
















Ten iterations of the algorithms are performed with a range of allowed dwell 
times to roughly illustrate the effect on localization error. The error is due to motion of 
the local node over the duration of the transmission, as described in Chapter IV. The 
trend in Figure 14 is similar to that expected when the vehicle speed is increased or 
processing time at the modems is increased. From ten iterations, it appears that neither 
the FP nor the DLC method is affected more than the other by this type of error. 

4. Data Sets with Nodes in a Straight Line 

As described in Chapter V, if there are only three nodes, the difference 
linearization algorithm must use all available data, even in cases with a bad sensor 
network geometry. If the nodes are in a straight line, the DLC method gives no solutions. 
The FP algorithm yields two ambiguous solutions for each waypoint. If neither solution 
has less error than the other, the algorithm arbitrarily chooses a solution. Half of the 
solutions will be correct and the other half will be mirror images of the correct solution, 
reflected about the axis on which the nodes are collinear, as shown in Figure 15. In this 
particular case, further filtering and comparison to previous positions can yield the 
correct position. For this thesis, the algorithms are limited to using only the current range 
information from the nodes. The previous vehicle position, speed and orientation are not 
used. Of note is the fact that the pairwise method yields solutions where the difference 
linearization method yields none. Figures showing the DL and FP method outputs for the 
same track, but with reference nodes in an approximately straight line, are included in 
Appendix F. 


33 



Actual and Calculated Vehicle Path for Simulated Data 



Figure 15. Nodes in a straight line will yield ambiguous solutions using the FP method. 

Incorporating information about the previous vehicle position, speed and 
orientation can resolve the ambiguity. 


Most schemes for node deployments, including that in [9], have reference nodes 
arranged in a straight line. In practice, this is not usually achievable. Nevertheless, 
deploying nodes in a straight line should be actively avoided if they are to be used for 
localization. 


5. Data Sets Varying the Number of Nodes 

Synthetic network geometries are generated with three to six remote nodes and 
the average localization error is evaluated for the three algorithms. The average 
localization error decreases as more remote nodes are made available for ranging, as 
shown in Figure 16. 
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Localization Error for 200 Iterations of Error-Free 
Synthetic Data, Varying the Number of Nodes 
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Figure 16. Average localization error for the three methods, varying the number of 

remote reference nodes. 


In overdetermined cases, the DLC method performs better than the FP method. If 
simulations are run for 3 remote nodes in the domain, the DL method and DLC method 
have the same output, as no data can be discarded in the case of a bad geometry, and the 
FP method performs best. When data from only good geometries were analyzed for 100 
iterations, the DLC method has an average localization error of 0.2473 meters, whereas 
the FP method has an average error of 0.3007 meters. For the randomly generated node 
positions, 24% of the geometries had bad condition numbers and were ignored in this 
calculation. The DLC method therefore performs better than the FP method for error-free 
data with any number of nodes, given a good sensor network geometry. With a bad 
geometry and only three remote reference nodes, the FP method performs better than the 
DLC method. 
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VIII. DESCRIPTION OF EXPERIMENT 


A. EXPERIMENTAL PLAN OVERVIEW 

An experiment was conducted in Del Monte Lake on the Naval Postgraduate 
School campus on 3-6 April and 16 April 2012. An aluminum boat is driven on the lake 
using a battery-powered trolling motor. The boat tows a sled behind it, fitted with a 
Seaweb modem and two GPS loggers. The towed modem, or local node, is controlled via 
a serial connection to a laptop computer on the boat. Commands to the local node are sent 
via a Procomm terminal on the laptop. Six Seaweb modems are deployed as remote 
nodes in the lake at the positions shown in Appendix A. As the boat transits the lake, the 
GPS loggers record fixes and the local node is commanded to issue pings to the remote 
nodes. The data packets returned in the echoes are recorded on a capture file for later 
analysis. A description of the equipment used and the experimental plan are found in 
Appendix G. 

B. SHALLOW LAKE ENVIRONMENT 

Del Monte Lake is a shallow fresh-water lake, approximately 150 meters by 150 
meters. During the experiment conducted in April 2012, fathometer measurements 
showed the depth to vary between 4 and 8 feet. The lake is fed by surface runoff water 
and has a silt bottom. A fountain in the center of the lake can be turned off to reduce 
ambient noise. A drain pump runs constantly to circulate water through the lake. 


<r 


Surface float 



Modem fixed 
to towed sled 


i 


Telesonar 

modem 


Mushroom anchor 


Average Depth 
(~2m) 


Waveguide Length (ISO meters) 


Air 

p = 1.21 kg/m 1 
c = 343 m/s 


Freshwater 
p * 998 kg/m 1 
c * 1481 m/t 


Muddy Bottom Layer 
p = 1600 kg/m* 
e * 1563 m/a 


Figure 17. Diagram of fixed node and mobile node towed by a surface vehicle. Sound 
speeds, densities and the lake dimensions are estimated for modeling sound 
propagation in Del Monte Lake. 
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The lake’s small volume, shallow depth and the drain pump make the lake 
susceptible to sound speed variability. The temperature of the water is affected by the 
amount of sunlight. The second week of tests occurred after a week of rain, which 
changed the lake’s salinity. The influx of sediment and drain water also influences the 
salinity. The constant motion of water and sediment into and out of the lake affects the 
bottom composition and bathymetry. 

For this analysis, the temperature and salinity are assumed to be spatially constant 
in the lake, and thus the sound speed is constant with range and depth. The zero-valued 
gradient leads to an infinite radius of curvature, or straight-line unrefracted propagation. 

For the short ranges in the lake, the transmission loss due to geometric spreading 
and attenuation is neglected. Noise can be a limiting factor. Noise sources in the lake 
experiments included the motor noise, flow noise across the transducer, flow noise across 
the body of the boat and towed sled, noise from the pump and fountain, and noise due to 
wind. The wind increased the lapping noise of water against the edge of the aluminum 
boat and sled. This was significant due to the shallow depth and proximity of the local 
transducer. When active, the fountain created broadband noise that interfered with the 
measurements. When collecting data, the fountain was turned off and the boat was 
operated at low speeds to minimize motor and flow noise. Data were taken on days when 
the local wind speeds were low and the surface was calm. 

Multipath is a significant problem in the shallow lake where there are many 
reflections from the bottom and surface. The MSFK modulation of the acoustic modems 
is designed to tolerate the impact of multipath. The shallow lake constrains the distance 
from the bottom and surface that the transducers can be suspended. 

For initial estimates of sound propagation, the lake is modeled as a waveguide 
with a depth of 2 meters and a length of 150 meters. The upper boundary is represented 
as a pressure release boundary. Density and sound speed values of 1.21 kg/m and 343 
m/s are used for air. The lower bound is also a pressure release boundary. Typical density 
and sound speed values for a muddy bottom are 1600 kg/m and 1563 m/s. For freshwater 
a density of 998 kg/m and a constant sound speed of 1481 m/s are assumed. The remote 
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source depth is set to 1 meter and is the approximate average position of transducers on 
the anchored modems. The receiver depth of 0.5 meters represents the local transducer 
mounted under the towed sled. 


BELLHOP- Del Monte profile 
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Figure 18. Bellhop sound propagation modeling software is used to produce eigenray and 
impulse response plots given a sound speed profde, medium characteristics, 
source and receiver locations, and dimensions of the operating space. From [21]. 


The eigenrays for sound propagation in the lake are modeled using the Bellhop 
acoustic propagation model. Bellhop is a numerical ray-tracing program run via Matlab 
that can output ray amplitudes, paths, travel times, and transmission loss, given an 
estimated sound speed profile for the lake environment. A program description and user 
manual is provided in [22], Figure 18, created by Pascal Gagnon with Bellhop, shows a 
graph of the sound speed profile, the eigenrays for a finite number of launch angles, and 
the impulse response at the receiver for those eigenrays. The ray plot and impulse 
response show the amount of multipath in the lake. Of significance is the similarity in 
arrival time and amplitude of the direct path and the path with one surface reflection. The 
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impulse response shows that the amount of time it takes before the multipath signals die 
out in the channel is about 1 millisecond. 

C. IN-AIR TESTING OF EQUIPMENT 

Before the modems were deployed in the lake, the acoustic communications were 
tested in air. The Teledyne Benthos modems have a number of settings that can be 
adjusted for a specific mode of operation. Settings are easily changed in air at short 
distances, when the acoustic transmissions are fully reliable. For this experiment, settings 
were chosen to minimize the time between transmissions. The in-air testing was useful 
for familiarization with the equipment, the Procomm coding software, and the Seaweb 
commands and data packets. 
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IX. PRACTICAL RANGING APPROACHES 


A. BROADCAST PING 

During in-air testing, it was determined that the broadcast ping cannot be used for 
ranging in the lake experiments. The broadcast ping command is designed for ocean 
applications in which the nodes are separated by long distances and propagation times are 
longer. The default settings for a broadcast ping cause the local node to wait 1 minute for 
up to 16 responses. This time can be shortened to approximately 30 seconds if waiting for 
only six nodes. 

This is a problem for the planned tracking experiment. The 30-second wait time 
limits the update frequency at which waypoints can be calculated. In the ocean, where the 
distance between nodes can be 1-2 kilometers, a waypoint every 30 seconds is ideal. In 
the lake, the largest distance between the fixed and mobile nodes is approximately 200 
meters. The boat travels approximately 2 miles per hour and can therefore cover about 30 
meters in 30 seconds. 30 meters is a significant portion of the network domain. Further, 
all ranges calculated during the 30-second wait time are used to plot one waypoint. This 
and non-stationary node geometry lead to significant localization errors. Finally, though 
the wait time for the local node is reduced to 30 seconds, the dwell times for the remote 
nodes cannot be changed by the user. Not all remote nodes return an echo due to dwell 
times greater than 30 seconds. 

The lowest possible dwell and wait times are desired to calculate ranges that 
approximately represent the same position, and to allow for a high update frequency of 
waypoint calculations. 

B. SEQUENTIAL PINGING 

An alternative to the broadcast ping command is to use the addressed ping 
command between two nodes. This method was used for the lake experiments. This 
function does not allow the local node to communicate with all remote nodes 
simultaneously. Each ping is addressed to a specific remote node. With six remote nodes, 
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six separate pings need to be sent out, versus one with a broadcast ping. The echoes 
return immediately, however, which ultimately reduces processing time. 

There are several constraints on the dwell and wait times, including 
environmental limits, signal processing limits, and operating limits of the equipment and 
coding software. The round-trip propagation time limits the wait time for the local node. 
In the lake, the round-trip times are on the order of milliseconds and are a negligible 
limitation. The processing time at local and remote nodes can be reduced by adjusting the 
modem settings. Communications can be impaired if the modems issue signals before the 
incoming multipath signals die out. A delay is automatically built into the processing 
time of the modems to avoid multipath interference. The lowest total processing time 
achieved was 4.2 seconds. For six nodes, echoes can be collected in 25.2 seconds, which 
is an improvement from the 30 seconds required for the broadcast ping. If an echo is not 
received, the local node waits for a set amount of time before continuing with the next 
ping. The lowest wait time achieved by adjusting modem settings was 10 seconds. This 
causes large gaps in the data due to one bad transmission. The 4.2 seconds between 
outgoing pings yields an error equivalent to the vehicle motion error described in Chapter 
IV. This method of data collection is not ideal, but it is more reliable than the broadcast 
ping for the lake experiments. 

C. LONG BASELINE MODE 

A Long Baseline (LBL) mode of operation for acoustic modems will resolve 
issues with the broadcast and sequential pinging modes. Short and Ultra Short Baseline 
systems also exist. The LBL systems, described in [8], are designed for ranges of 50 
meters to 2 kilometers. They are used commercially in the same manner as the Seaweb 
sensor networks. The systems use signal propagation times to calculate ranges and locate 
an object. The difference between the two lies in the form of communications. The 
baseline systems use transponders that receive a signal in one band of frequencies and 
reply in a different band. The proposed adaptation of the LBL systems for Seaweb 
networks would have each remote node respond to a broadcast ping using a unique set of 
frequencies. This would allow the echoes to return to the local node simultaneously and 
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would minimize errors due to vehicle motion over the duration of the transmissions. It 
would also permit a higher update frequency for ranging calculations. 
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X. FURTHER EXPERIMENTAL CONSIDERATIONS 


A. SPEED OF SOUND CALCULATION 

While conducting the experiment on Del Monte Lake, ranging commands were 
sent between the stationary, remote nodes to estimate the speed of sound in the lake. GPS 
positions were recorded when the anchored modems were deployed and recovered. 
Assuming the GPS coordinates are the exact locations of the modems, we can calculate 
the range between remote nodes and then back-calculate the sound speed using measured 
delay times. Formulae for calculating the distance between GPS points can be found 
online. The Haversine Formula used by [23] calculates the great-circle distance between 
two points, assuming a spherical earth with a radius of 6,371 kilometers. 

Coordinate-to-distance conversions depend on the particular latitude and 
longitude of the points. Choosing the coordinate position for Node #21 in the middle of 
the lake, located at 36°35.953’N and 121°52.203’W, one can find a conversion factor by 
calculating the distance between this point, and a point that varies by a given number of 
degrees north and west. The conversion factors were found as shown in Table 3, using 
the fonnula from [23]. 



Deg/Min N 

Deg/Min W 

Degrees N 

Degrees W 

Range (m) 

Node 21 

36°35.953’ 

121°52.203’ 

36.59922 

121.87005 


Point 2 



36.59923 

121.87005 


A° 



0.00001 

0.0 

1.112 



Deg/Min N 

Deg/Min W 

Degrees N 

Degrees W 

Range (m) 

Node 21 

36°35.953’ 

121°52.203’ 

36.59922 

121.87005 


Point 2 



36.59922 

121.87006 


A° 



0.0 

0.00001 

0.8927 


Table 3. A conversion factor between degrees and distance, specific to the operating area, 
is found by choosing two points separated by 0.00001 degrees in either direction. 


The conversion factors are used to calculate the ranges between each combination 
of nodes given their GPS coordinates. The ranges shown are also corrected for the depth 
of the nodes. Reversing the calculations in Chapter IV, the node-to-node ranges are 
converted from two dimensions to three. 
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The transmission delay times were collected between nodes on 16 April 2012. 
The sound speed is found by dividing the GPS-derived range by half the round-trip delay 
time. Table 4 shows the ranges, delay times and sound speeds calculated for nine node 
combinations. 


Node # - 
Node# 

Ranges for 
16 APR (m) 

Round-trip Delay 
for 16 APR (s) 

Sound 
Speed (m/s) 

20-21 

140.21 

0.1988 

1410.56 

20-23 

127.56 

0.1835 

1390.31 

20-25 

164.07 

0.2281 

1438.62 

21-22 

123.45 

0.1725 

1431.32 

21-23 

120.13 

0.1645 

1460.60 

22-24 

103.08 

0.1438 

1433.62 

22-25 

123.24 

0.1691 

1457.65 

23-24 

134.44 

0.185 

1453.39 

24-25 

114.80 

0.1552 

1479.42 

Average 



1439.50 


Table 4. The speed of sound in the lake is found using ranges between anchored modems. 

This method yields a large range of solutions due to errors in the measurements. 


The average calculated sound speed is 1439.5 m/s with the individual sound 
speeds varying over a range of 89.11 m/s. The average speed is far from the 1481 m/s 
sound speed estimate for freshwater, but within the range of typical sound speed values 
for freshwater. 

B. GPS ERROR 

Table 5 contains the GPS coordinates recorded on 3-4 April 2012, when the 
nodes were deployed, and those from 16 April 2012, recorded during recovery. 


Date 


Node 20 

Node 21 

Node 22 

Node 23 

Node 24 

Node 25 

3-4 APR 

N 

36°36.009’ 

36°35.953’ 

36°36.018’ 

36°36.015’ 

36°35.962’ 

36°35.966’ 


W 

121°52.270’ 

121°52.203’ 

121°52.228’ 

121°52.182’ 

121°52.246’ 

121°52.171’ 

16-Apr 

N 

36°36.005’ 

36°35.951’ 

36°36.015’ 

36°36.014’ 

36°35.962’ 

36°35.966’ 


W 

121°52.268’ 

121°52.202’ 

121°52.225’ 

121°52.183’ 

121°52.246’ 

121°52.169’ 


Table 5. GPS coordinates for the anchored modems when deployed and recovered. 
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Table 6 compares the ranges between the nine sets of nodes for the deployment and 
recovery dates. A comparison of the initial and final GPS locations shows that the nodes 
remained approximately stationary throughout the testing, however differences in the 
GPS measurements correspond to range differences of up to 5.13 meters between the 
same two nodes. From visual inspection, it was clear that the anchored modems were not 
displaced by this amount. The range difference can be attributed to error in the GPS 
measurements. 


Node # - 
Node# 

Ranges for 
3-4 APR (m) 

Ranges for 
16 APR (m) 

|AR| (m) 

20-21 

143.91 

140.21 

3.70 

20-23 

131.40 

127.56 

3.84 

20-25 

167.47 

164.07 

3.40 

21-22 

126.08 

123.45 

2.63 

21-23 

119.08 

120.13 

1.06 

22-24 

107.19 

103.08 

4.11 

22-25 

128.37 

123.24 

5.13 

23-24 

136.80 

134.44 

2.37 

24-25 

111.83 

114.80 

2.97 


Table 6. Ranges between anchored modems are calculated using their GPS coordinates. 


Using the calculated average sound speed reduces the localization error for both 
algorithms compared to the error found with the typical freshwater sound speed of 1481 
m/s. The error in the GPS readings significantly affects the calculated sound speed. For 
this reason, a secondary method for calculating sound speed is desired. A common 
practice is to use a conductivity-temperature-depth (CTD) probe to directly measure the 
water properties that determine sound speed. 
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GPS Tracks 



Figure 19. Plots of the two simultaneously logged GPS tracks shows that significant error 
can exist in the reference track. Tracks with good agreement were chosen for data 

analysis. 

The modem positions were recorded with a Garmin GPSmap76CS logger. The 
two G-Log 760 GPS loggers used to record the vehicle track also had errors, as seen by 
plotting their output. The GPS loggers were placed approximately half a meter apart, but 
recorded positions up to 20 meters apart. For this reason, the tracks used to analyze the 
error in the localization algorithms were limited to those with good agreement between 
the two GPS loggers. 
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XI. ALGORITHM IMPLEMENTATION FOR EXPERIMENTAL 

DATA 

A. PROGRAM DESCRIPTION 

1. Initial Data Processing 

GPS fixes are retrieved from the loggers using the accompanying software and 
saved as Microsoft Excel files. The displayed data from node-to-node communications, 
shown in Figure 3, are saved by Procomm in a text-formatted capture file. A Matlab 
script reads in the GPS data. For consistency, the GPS coordinates of the nodes recorded 
on 16 April 2012 are used when analyzing all tracks. For ease of plotting and 
computation, the latitude and longitude fixes of the GPS-logged tracks, and the locations 
of the remote nodes are converted to x and y coordinates with units of meters and with the 
origin defined to be a set position southwest of the operating area in Del Monte Fake. 
The x axis runs west to east and the y axis runs south to north. Again, the conversion 
factors between GPS coordinates and meters were found as shown in Table 3, using the 
formula from [23]. The text files with range data are sorted through with a second Matlab 
script to extract the required infonnation. 

After GPS data and data from Procomm capture files are converted to Matlab data 
files, they are read into a Matlab script specific to each algorithm, where the roundtrip 
propagation times are converted to ranges. The depths of the transducers are then used to 
project the ranges onto a two-dimensional plane. 

The track is manually truncated to only include lengths of track for which there is 
a sufficient amount of recorded node-to-node ranges. For the experiment, the individual 
remote nodes were pinged using a programmable hotkey via Procomm. An infinite loop 
script was less efficient due to irregular transmission times and errors. Each modem can 
only store a given number of bits in a temporary cache. If the local node is commanded 
via a hotkey to send a ping before the incoming ping is processed, the cache will fill and 
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the modem either pauses or automatically reboots. The sections of sparse data shown in 
Figure 20 occurred when transiting around the fountain or when sequential pings were 
sent too quickly. 


GPS Tracks and Calculated Track 
for 6 Ranges per Waypoint with FP Method 



Figure 20. The FP algorithm output is shown for data taken on 16 April 2012. Sections of 
sparse data occurred when transiting around the fountain or when sequential pings 

were sent too quickly. 


Both the GPS data and Seaweb data are recorded with Universal Time 
Coordinated (UTC) time. It is assumed that the times from both GPS loggers time and the 
local node are correct. This cannot be exactly verified using the equipment. Data from the 
local node are passed to the Procomm terminal that has its own clock, based on the 
computer clock. The time it takes to pass this information will skew any calibration. The 
GPS loggers used do not allow data to be viewed as they are recorded. Imprecise methods 
of checking the time agreement show that the times do not differ more than a few 
minutes, thus the assumption is maintained for initial algorithm testing. The GPS track is 
truncated to correspond to the Seaweb data times and the times are shifted to begin at 
zero. 
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2. Sorting Data 

The necessary use of sequential pinging as the vehicle moves introduces 
significant error when using trilateration. An example of this error is shown in Figure 21. 


GPS Tracks and Calculated Track 
for 6 Ranges per Waypoint with Pairwise Method 


GPS Tracks and Calculated Track 
for 6 Ranges per Waypoint with Pairwise Method 



# Points along 
vehicle track at 
which echoes 
arrive, labeled 
with the node 
number 


-—GPS Path 

A Position of submerged nodes 
X Calculated avg position 
-Range circles 



E-W Distance (m) 


E-W Distance (m) 


Figure 21. The first image shows range circles from six sequential pings, the vehicle 
track, and the average position calculated with the filtered pairwise algorithm with 
the six returns. The second image shows the same at a higher resolution. 


The figure shows that the returns are spaced almost evenly along the track. They 
are each separated by approximately 4.2 seconds. The even spacing indicates that little 
error exists in the individual node-to-node range data. The average position is calculated 
using the range circle crossings, most of which are shown. One can see that the range 
circles from Node 1 and Node 4 do not intersect, as well as for Node 1 and Node 3. For 
the first set, the approximate solution described in Chapter IV is a good approximation, 
whereas it is worse for the second node set. 

To minimize error due to vehicle motion, and maximize the amount of waypoints 
calculated, the range data are broken up in overlapping sections of lengths 3 through 6, as 
shown in Figure 22. 
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Figure 22. The schematic shows successful echo returns from nodes 1-5. The return from 

node 6 was unsuccessful. Two echoes from node 3 were recorded in a row. These 
errors limit the division of data for analysis. 

Data sections with three ranges are closer in time. Data sections with six ranges 
will be more spaced out, but will allow for more averages. If an echo is not received from 
a node, the returns before and after the data gap will not be used to calculate a waypoint. 
This is due to the fact that the local node waits for a return for approximately ten seconds 
before continuing with the next ping, yielding a total of about 15 seconds between ranges. 
On a few occasions, pings to the same node were sent in sequence. The two ranges are 
unusable, as they will yield range circles that are an equal distance apart in all directions. 
These errors limit the sections of six ranges more so than the sections of three ranges. 
The effect is evident in the plotted results included in Appendix B. 

3, Determining Average Localization Error 

After each Matlab script sections the data, the filtered pairwise (FP) and the 
difference linearization with a condition check (DLC) localization algorithms are 
implemented to calculate the average position for each waypoint. The localization error is 
found by computing the distance in meters from each calculated waypoint to the 
corresponding waypoint of the GPS track. All distances are summed and then divided by 
the number of waypoints used to find the average localization error for the track. For 
consistency, when calculating the localization error for each algorithm, the same number 
of waypoints is used for the same length of track. 

Difficulties arise when trying to find the true position and the corresponding 
calculated waypoint due to GPS logger error, the inability to synchronize GPS logger 
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time with Seaweb node clocks and Procomm time, the resolution of GPS tracks, and 
interpolation in Matlab, and the fact that the return echoes used to calculate one waypoint 
occur at different times. 

The time chosen to represent each waypoint is the median time that the echoes 
return. The Matlab timeseries function is used to interpolate GPS track positions at these 
median times. The positions from the two GPS tracks are averaged at each time for 
comparison with the algorithm output. The GPS loggers are set to take one fix every 
second. The interpolation error in Matlab is minimal. 

B. EMPIRICAL CALIBRATIONS 

Once the algorithms output an average localization error for the track, a time 
offset is added to the Seaweb echo times to align them with the GPS UTC time, and the 
average localization error for the track is minimized. 


Time Offset Calibration for 6 APR 2012 Track 


Time Offset Calibration for 16 APR 2012 Track 




Figure 23. Calibration for the time offset between Seaweb clocks and GPS logger clocks. 


Testing the algorithm over a range of time offsets, it is found for the April 6 data 
that the Seaweb clocks are approximately 57.5 seconds ahead of the GPS times, and 61 
seconds ahead for the April 16 data, as shown in Figure 23. For the calibration shown, the 
algorithms are run with sets of three ranges per waypoint. Further examination shows that 
the error is minimized for any sized set of ranges for these time offsets. 
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Figure 24. Determination of the sound speed that minimizes localization error. 


After the time offset is applied to the data, a similar calibration is done to find the 
sound speed that yields the least average error. For the calibration shown as shown in 
Figure 24, the algorithms are run with sets of six ranges per waypoint, as this setting 
yields the minimal localization error at any sound speed. The sound speeds are expected 
to differ for the two days. For the April 6 data, a sound speed of approximately 1447.75 
m/s minimizes the localization error. For the April 16 data a sound speed of 
approximately 1426 m/s minimizes the error. 

C. RESULTS DISCUSSION 

The plotted results from the algorithms are included in Appendix B and 
summarized in Table 7. Both algorithms performed best when using combinations of six 
ranges to find each waypoint. Though the vehicle positions are calculated using echo 
returns over a range of approximately 25 seconds, versus about 13 seconds for sets of 
three ranges, the amount of averaging allowed with more range information has a greater 
effect on algorithm performance. 


Average Localization Error in Meters for Each Method 


Data from 6 APR 2012 

Data from 16 APR 2012 

Ranges Used 

DLC Method 

FP Method 

DLC Method 

FP Method 

3 

4.99 

4.12 

5.30 

4.34 

4 

4.64 

3.63 

4.93 

3.95 

5 

4.14 

3.32 

4.52 

3.75 

6 

3.62 

3.04 

4.20 

3.65 


Table 7. Average localization error for the FP and DLC methods for experimental data. 
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The plotted tracks show erratic discontinuities that appear to be semi-periodic in 
localized regions. These discontinuities are not due to errors in the data, but are results of 
sequential pinging. The semi-periodicity is a result of using overlapping sections of data 
for nearby waypoints. The synthetic simulations were run assuming an LBL-type mode 
of operation for which a single ping communicates with all remote nodes. The remote 
nodes reply with no delay and can arrive at the local node simultaneously. If the synthetic 
data are made to mimic the sequential pinging mode of communications, plots similar to 
that in Figure 25 are produced. 


Synthetic Track and Calculated Track 
with FP Method and Sequential Pinging 
for 5 Ranges per Waypoint 



Figure 25. Example of discontinuities over a synthetic track due to sequential pinging. 

The figure shows the same discontinuities over the randomly generated synthetic 
track as seen in the experimental data. The simulation was run for error-free synthetic 
data. There were no delays added at the remote nodes. Pings were sent once every 3 
seconds. Even with the high ping frequency and no error, the pairwise method yields an 
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average error of about 3 meters. The average localization error of the FP and DLC 
methods with experimental data cannot be expected to be better than 3 meters. 

The FP algorithm outperformed the DLC algorithm regardless of the number of 
reference nodes. This can be explained by the FP method’s robustness with regards to 
range error. In most simulations, the DLC algorithm has a lower average localization 
error, though by a small margin. Figure 13 shows that when error in the individual range 
measurements is present, the FP method performs better than the DLC on average. When 
pinging remote nodes sequentially, the displacement of the vehicle over the period of 
time can translate to individual range measurement errors. 

D. CONDITION NUMBERS REVISITED 

In looking at the eigenvalues and condition numbers of the matrices for the 
deployed nodes, it was discovered that each algorithm had different condition numbers 
using the same node locations. This led to the realization that for every three nodes, there 
are in fact three condition numbers possible. 

p f(x 2 -Xj) 

U X 3-X 2 ) (P 3 T 2 )y 

Q J(x 3 —Xi) (w-Tj)" 

^(x 3 -x 2 ) (y 3 -y 2 ), 

R f(x 3 —x,) (t 3 -Ti)" 
l(Xl-X 2 ) (Ti -p 2 )y 

The A matrix from Equation (23) can be expressed as one of the three matrices 
shown, each of which has a different condition number. Of the modems deployed in the 
lake, Nodes 20, 22 and 23, pictured in Appendix A, are approximately in a straight line, 
and their A matrices yield the largest condition numbers. The three condition values are 
16.16, 15.74 and 6.38. The condition numbers for all A matrices for all node 
combinations are included in Appendix C. 

The DLC algorithm can be edited to compute the three condition numbers and 
find the minimum of the three. The range equations will then be summed in the way that 
corresponds to the lowest number. Editing the DLC code in this way approximately 
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doubles the computation time. It also does not guarantee that the error will be reduced. 
The calculated condition of a matrix is an upper limit, and as seen in Equation (30), the 
fractional error in a solution is less than or equal to the condition number multiplied by 
the fractional error in the measurement. In Table 8, the average localization errors are 
compared for the FP method, the DLC method using the P matrix, and the DLC method 
using the minimum condition value. 


Average Localization Error in Meters for Each Method 


Data from 6 APR 2012 

Data from 16 APR 2012 

Ranges 

DLC Method 

with P Matrix 

DLC Method 
Min Condition # 

FP 

DLC Method 

with P Matrix 

DLC Method 

Min Condition # 

FP 

3 

4.99 

4.99 

4.12 

5.30 

5.30 

4.34 

4 

4.64 

4.16 

3.63 

4.93 

4.81 

3.95 

5 

4.14 

4.37 

3.32 

4.52 

4.78 

3.75 

6 

3.62 

3.63 

3.04 

4.20 

4.41 

3.65 


Table 8. Average localization error for experimental data using the FP method, the DLC 
method using one condition number, and the DLC method that uses the minimum 
of all three possible condition values of the inverted A matrix. 


The table shows that using the matrices with the lowest condition numbers does 
not necessarily reduce the overall localization error for the experimental data. When 
using three ranges, the average error is the same to 14 decimal places, for example. The 
Node 20, 22, and 23 combination conveniently never occurs when using three ranges due 
to sequential pinging. The error was only reduced when using combinations of four 
ranges. The plotted results are shown in Appendix D. The FP method still has the lowest 
error in all cases. 

Finally, the only area in which the FP method performed better than the DLC 
method in simulations was for synthetic data with error introduced to each individual 
range measurement. To verily that the FP algorithm performs better than the DLC 
method using the smallest condition numbers, 50 iterations are repeated with the 
corrected DLC method. 
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Figure 26. Average localization error as it increases with the amount of error introduced 
to each individual range measurement. The FP method is compared to the DLC 
method that uses the minimum matrix condition number. 


In Figure 26 it is seen that the localization error using the FP method still grows 
at a slower rate than for the DLC method. Thus, making this final correction to the DLC 
method does not guarantee that the localization error will be reduced. 
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XII. PRACTICAL IMPLEMENTATION OF LOCALIZATION 

ALGORITHMS 

Though the simulations are run for ideal data and the experimental data are taken 
in a controlled environment, the development of the filtered pairwise (FP) and the 
difference linearization with a condition check (DLC) localization algorithms is done 
with practical implementation in mind. In evaluating localization algorithms, important 
performance metrics are the computation time, the robustness to error in the ranges, and 
the operational limits and cost of the existing equipment. 

A. DEPLOYMENT CONSIDERATIONS 

Modems deployed in the ocean will typically be as sparsely distributed as 
possible, so as to cover the greatest area with minimal cost. Trilateration requires the 
communication area of three neighboring nodes to overlap. Only three nodes are 
available for ranging in most cases, versus the six nodes used in this thesis, thus it is 
critical that those three nodes are not deployed in a geometry that results in bad matrix 
condition numbers. 

B. REAL-TIME IMPLEMENTATION CAPABILITY 

The computation time of each algorithm is an important consideration for real¬ 
time implementation. The computation times for the experimental data are summarized in 
Table 11. The times shown are for the computation of the entire track. Depending on the 
length of the track, the algorithm used, and the number of ranges per waypoint, each 
algorithm performs a different number of calculations. Table 9 shows the number of 
waypoints calculated for each track, varying with the number of range data used to find 
each waypoint. The data from the April 6 th track has a greater decrease in calculated 
waypoints as the ranges per waypoint increase, than the data from the April 16 th track. 
This can be attributed to the sparseness of data. 
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Number of Waypoints Calculated for Both FP and DL Methods 

Ranges per Waypoint 

Data from 6 APR 2012 

Data from 16 APR 2012 

3 

239 

364 

4 

214 

352 

5 

189 

341 

6 

162 

330 


Table 9. The number of waypoints calculated for each track depends on the length of the 
track and number of range data used to find each waypoint. 


Table 10 shows the total number of combinations used to calculate the waypoints 
for each track and for each localization algorithm, as the number of ranges per waypoint 
increase. The numbers of waypoints from Table 9 are multiplied by the number of 
combinations of range circles used to find each solution. The multiplicative factors are 
calculated with the choose k' formula in Equation (32). The DL and DLC methods use 
combinations of three range data, whereas the FP method uses combinations of two range 
data to find a solution. This difference, combined with the number of ranges used to 
calculate each waypoint, is the greatest contribution to algorithm computation time. 


Number of Range Data Combinations Used 
for Calculation of Waypoints over Entire Track 


Data from 6 APR 2012 

Data from 16 APR 2012 

Ranges per wp 

DL Method 

FP Method 

DL Method 

FP Method 

3 

239x 1 = 239 

239x3 = 717 

364 x 1 = 364 

364x3 = 1092 

4 

214x4 = 856 

214x6 = 1284 

352x4= 1408 

352x6 = 2112 

5 

189 x 10 = 1890 

189 x 10 = 1890 

341 x 10 = 3410 

341x10 = 3410 

6 

162 x20 = 3240 

162x 15 = 2430 

330x20 = 6600 

330x15 = 4950 


Table 10. The DL method uses combinations of three range data, whereas the FP method 
uses combinations of two range data to find a solution. The difference has a large 
effect on algorithm computation time. 


If five range data are used per waypoint calculation, both algorithms have the 
same amount of combinations to compute. Thus, the computation times compared for this 
case reflect the complexity of each algorithm’s calculations. Table 11 shows that the FP 
method has the lowest computation time per combination. The multiplicative factors in 
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Table 10 are also the number of averages that each method can use to calculate a 
waypoint, assuming no data is discarded from bad geometries. 


Computation Time in Seconds for Each Method 


Data from 6 APR 2012 

Data from 16 APR 2012 

Ranges 
per wp 

DLC Method 

with P Matrix 

DLC Method 

Min Condition # 

FP 

DLC Method 

with P Matrix 

DLC Method 

Min Condition # 

FP 

3 

0.033 

0.053 

0.066 

0.046 

0.079 

0.085 

4 

0.077 

0.139 

0.083 

0.121 

0.223 

0.113 

5 

0.150 

0.268 

0.086 

0.260 

0.475 

0.138 

6 

0.244 

0.448 

0.103 

0.487 

0.889 

0.186 


Table 11. Comparison of computation times for the FP and DLC methods using 

experimental data, depending on the number of ranges used to calculate each 

waypoint. 


Table 11 also shows that the DLC method that minimizes the matrix condition 
numbers has a computation time almost twice that of the DLC method using only the P 
matrix. The computation times in the table do not include the time needed to preprocess 
or post-process the data. Again, in practical situations, only three nodes are typically 
available for ranging. For real-time implementation of the localization algorithms, the 
trade-off between computation time and the amount of averaging enabled by each 
algorithm should be considered. 

C. CRITICAL ASSUMPTIONS 

Assumptions about the controlled lake environment will not hold in an ocean 
environment. In particular, the constant sound speed and two-dimensionality 
approximations are critical assumptions that become less accurate as the range and depth 
of the operational domain increase. Additional complexity of the environmental model is 
required, but is hampered by computation time considerations and spatially and 
temporally changing parameters. For this reason, a robust localization algorithm is 
required, versus one that requires near-ideal data to converge on a solution. 
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D. AVAILABILITY OF ADDITIONAL INFORMATION 


Finally, when operating a sensor node on a moving vehicle, it is likely that the 
vehicle’s heading and speed will be known. One might also kn ow the depth of the vehicle 
and the local sound speed. This additional infonnation can aid in localization by reducing 
or constraining the degrees of freedom. 

1. Vehicle Speed and Heading 

A vehicle’s speed and heading are, of course, helpful in determining its future 
location. The plot in Figure 27 was created by iteratively solving a system of nonlinear 
range circle and distance interval equations, assuming a constant heading and a constant 
time interval between pings. 


Waypoints Calculated for Straight Line Path 
with Known Vehicle Speed 



E-W Distance (m) 


Figure 27. A vehicle’s speed can be used to solve for its future positions. 

Initial testing of this method, inspired by Ensign Michael Moberg, showed that for 
ideal synthetic data, the iterative nature of the function solver yielded inaccurate solutions 
in about 50% of the simulations. When the function solver converged on a correct 
solution, the results were highly accurate. This iterative approach was computationally 
inefficient and was not further investigated, but it serves as a simple example of 
incorporating additional information into the localization algorithm. 
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2 . 


Previous Positions 


The previous position of a vehicle is valuable information, and can be used to 
filter out ambiguous solutions. In many cases, the incorrect ambiguous solution from 
intersecting range circles is more distant from the previous vehicle position than is 
physically possible, given the vehicle’s speed. Further, if the pings are sent at a sufficient 
frequency or the track of the vehicle is unchanging, the previous two positions of the 
vehicle can be used to estimate the vehicle heading. 
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XIII. CONCLUSIONS AND RECOMMENDATIONS 
FOR FUTURE WORK 

A. THESIS CONCLUSIONS 

Simulations of a mobile sensor node moving in the domain of a fixed sensor 
network were created. Synthetic node-to-node range data were generated to test existing 
range-based localization algorithms. 

Simulations showed that difficulties arise when using the previously developed 
difference linearization (DL) algorithm due to matrix inversion. When solving linear 
matrix equations, the condition number of a matrix characterizes the solution’s sensitivity 
to error in the equation inputs. The inverted matrix in the DL algorithm includes only the 
fixed positions of remote nodes. Thus, the condition number can be used to evaluate good 
network geometries. For bad network geometries, the localization error is large for the 
DL method. A correction to the difference linearization method is proposed that checks 
the condition of each inverted matrix and discards solutions calculated with ill- 
conditioned matrices. Difference linearization with a condition check (DLC) greatly 
improves the perfonnance of the DL method in over-detennined situations. 

The filtered pairwise (FP) algorithm developed in this thesis is a good alternative 
to the DL and DLC methods previously developed. The FP method does not require 
matrix inversion. Comparisons of the two localization algorithms in simulation showed 
that for good network geometries, the DLC method yielded lower average localization 
error with error-free data, when varying the number of remote nodes, and when dwell 
times are added at the remote nodes. The FP method is shown to be more robust when 
error is introduced to the individual range measurements. 

In practice, only three nodes will typically be available for ranging. If the three 
nodes are arranged in a geometry that yields ill-conditioned matrices, the FP method 
performs better than the DL and DLC methods. 

The sensor network deployed in Del Monte Lake on the NPS campus in April 
2012 was used to collect node-to-node ranges between the submerged, anchored remote 
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modems and a local modem towed behind a boat operating on the surface. Experimental 
data were collected and analyzed using the FP and the DLC algorithms. GPS fixes were 
recorded simultaneously to provide a reference track for analysis. 

During in-air testing of equipment, it was determined that the broadcast ping 
command was not suited for the lake experiments. The experimental data were taken by 
sequentially pinging individual remote nodes with the mobile local node. A simulation of 
sequential pinging with error-free synthetic data shows that the average localization error 
with the FP method has a lower limit of approximately 3 meters. 

In experimentation, the FP method calculated the track of the mobile node with 
more accuracy than the DEC method and had a lower computation time. 

The superior performance of the FP method over the DEC method with only three 
reference nodes available and the algorithm’s robustness make the simple FP method an 
attractive option for localization. The FP algorithm can be further improved by 
incorporating more information about the mobile vehicle, such as its speed, heading and 
depth. 

B. FUTURE WORK 

1. Kalman Filters 

Kalman filters are commonly used for tracking and navigation of autonomous 
vehicles with uncertain or noisy measurements. A Kalman filter algorithm is used to 
iteratively estimate the state of a linear system that is subject to Gaussian noise. The filter 
estimates a parameter vector based on previous measurements of the system, and 
minimizes the mean squared error of the estimate. New information is fed into the 
algorithm as it is received. A forgetting factor is used to truncate the data from previous 
measurements that are fed into the Kalman filter. The basic Kalman filter, forgetting 
factor and sample coded applications are described in [24], 
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time k-1 


step through time 


Figure 28. Schematic of a Kalman filter procedure. After [25]. 


An Extended Kalman Filter (EKF) is used for non-linear systems. The application 
of an EKF to localization is described in [25]. For a two-dimensional scenario, the 
system, or behavior of a vehicle, can be described with a linear matrix equation derived 
from basic kinematic equations. Shareef and Zhu [25] step through the procedure of an 
EKF algorithm for a trilateration problem with the range equations expressed in Equation 
(20). An EKF is used to make a position estimate, a position-velocity estimate, and a 
position-velocity-acceleration estimate depending on the kn own information. The benefit 
of using additional information depends on the particular scenario. The formulation in 
[25] can be extended to create a position-velocity-bearing model of an EKF. Ideally, the 
depth of a vehicle can be used for a three-dimensional model of an EKF. An Extended 
Kalman Filter can be tested using the experimental data already collected and no 
additional information. Further simulations and experimentation can investigate the 
benefit of adding speed and acceleration information, as discussed in [25], as well as 
other parameters. 

2. Experimental Controls 

When conducting future experiments in Del Monte Lake, further experimental 
controls are possible. The uncertainty in the sound speed can be resolved by taking 
conductivity-temperature-depth (CTD) probe measurements to determine the sound- 
speed-dependent medium characteristics. More accurate GPS loggers can be used, or 
their agreement can be checked prior to conducting lake tests. The GPS logger agreement 
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depends in part on the satellites that the loggers choose to use as references. The satellites 
can be explicitly selected with the accompanying software. 

A thorough survey of the lake, specifically its bathymetry and the bottom type, 
would aide in modeling the sound propagation. The remote nodes could also be more 
accurately located. The weights should be anchored and monitored to ensure their 
locations do not drift over time. The tethers can be better rigged to limit the watch circles 
of the underwater nodes. Recordings of the acoustic signals received at the local node can 
be studied to detennine the effect of multipath on the accuracy of the range data. Finally, 
a statistical analysis of the accuracy of range data can be conducted by ranging between 
stationary anchored nodes over a period of time. 

The environmental characterization and added experimental controls will enable 
better quantization of the localization error based on range measurements. 

C. APPLICATIONS 

The node-to-node range data for Seaweb acoustic sensor networks are collected as 
a by-product of normal acoustic communications. Accurate localization and tracking of 
mobile sensor nodes can be the basis for deployable underwater ranges in the future. 

When submerged, UUV’s and submarines must use inertial navigation to estimate 
their position. Fix expansion occurs over time, requiring the vehicles to surface and 
obtain a GPS fix. A deployable range of nodes can serve as an aide to submerged 
navigation. Submerged or surface vehicles fixed with a local node can use node-to-node 
range data to establish a position estimate in a GPS-denied environment. 

Localization algorithms developed for mobile nodes can also be implemented in a 
network to localize the stationary nodes. In deep water it can be difficult to accurately 
measure the location of nodes. Self-localization in a network is desired to initialize and 
update the positions of remote network nodes. 

The development of further practical applications for Seaweb sensor networks 
will benefit from accurate, robust, and real-time range-based localization algorithms. 
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APPENDIX A: DEL MONTE LAKE AND MODEM LOCATIONS 



The Google Earth image shows the approximate locations of the deployed modems in 
Del Monte Lake on the NPS campus. The GPS coordinates from the Garmin 
GPSmap76CS logger were entered into Google Earth. The pictured location of node 22 
indicates an error either in the logged GPS locations or in plotting the points via Google 
Earth. The node numbers refer to the node ID numbers used for node-to-node 
communications. Nodes 20-25 are referred to throughout this thesis as nodes 1-6. The ID 
number for local node, fixed to the towed sled, is 26. 
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APPENDIX B: PLOTTED EXPERIMENTAL RESULTS 


Contact professor J.A. Rice atjarice@nps.edu for a copy of the experimental data and the 
Matlab codes used to implement the localization algorithms and create these plots. 


GPS Tracks and Calculated Track 
for 3 Ranges per Waypoint with DLC Method 



GPS Tracks and Calculated Track 
for 3 Ranges per Waypoint with FP Method 



71 























N-S Distance (m) N-S Distance (m) 


GPS Tracks and Calculated Track 
for 4 Ranges per Waypoint with DLC Method 



GPS Tracks and Calculated Track 
for 4 Ranges per Waypoint with FP Method 
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N-S Distance (m) N-S Distance (m) 


GPS Tracks and Calculated Track 
for 5 Ranges per Waypoint with DLC Method 



GPS Tracks and Calculated Track 
for 5 Ranges per Waypoint with FP Method 
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N-S Distance (m) N-S Distance (m) 


GPS Tracks and Calculated Track 
for 6 Ranges per Waypoint with DLC Method 



GPS Tracks and Calculated Track 
for 6 Ranges per Waypoint with FP Method 
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N-S Distance (m) N-S Distance (m) 


GPS Tracks and Calculated Track 
for 3 Ranges per Waypoint with DLC Method 



GPS Tracks and Calculated Track 
for 3 Ranges per Waypoint with FP Method 
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N-S Distance (m) N-S Distance (m) 


GPS Tracks and Calculated Track 
for 4 Ranges per Waypoint with DLC Method 



GPS Tracks and Calculated Track 
for 4 Ranges per Waypoint with FP Method 
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N-S Distance (m) N-S Distance (m) 


GPS Tracks and Calculated Track 
for 5 Ranges per Waypoint with DLC Method 



GPS Tracks and Calculated Track 
for 5 Ranges per Waypoint with FP Method 
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N-S Distance (m) N-S Distance (m) 


GPS Tracks and Calculated Track 
for 6 Ranges per Waypoint with DLC Method 



GPS Tracks and Calculated Track 
for 6 Ranges per Waypoint with FP Method 
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APPENDIX C: CONDITION NUMBERS OF MATRICES 


Condition numbers 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Node Combos 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

21 

21 

21 

21 

22 

22 

22 

23 

23 

24 

22 

23 

24 

25 

23 

24 

25 

24 

25 

25 

for P matrix 

3.99 

1.84 

5.16 

2.60 

6.38 

2.18 

2.22 

2.89 

1.42 

1.56 

for Q matrix 

2.54 

2.02 

5.78 

5.93 

16.16 

1.32 

4.01 

1.61 

3.43 

3.37 

for R matrix 

1.85 

1.47 

2.21 

3.67 

15.74 

2.81 

5.55 

1.86 

2.67 

4.03 


Condition numbers 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

Node Combos 

21 

21 

21 

21 

21 

21 

22 

22 

22 

23 

22 

22 

22 

23 

23 

24 

23 

23 

24 

24 

23 

24 

25 

24 

25 

25 

24 

25 

25 

25 

for P matrix 

2.14 

3.36 

4.26 

3.68 

4.38 

6.18 

3.24 

1.59 

1.49 

2.64 

for Q matrix 

3.77 

2.41 

2.28 

1.76 

3.25 

2.37 

1.79 

3.14 

1.78 

2.07 

for R matrix 

1.99 

1.50 

2.27 

2.34 

1.76 

5.62 

4.40 

4.01 

2.10 

1.33 


The P, Q and R matrices are defined as: 

p r(* 2 -*i) (T 2 -.>if 
^(x 3 -x 2 ) (y 3 ~y 2 )^ 

j o 3 -*,) (T3-T1)" 
^(x 3 -x 2 ) (y 3 -y 2 ) y 

R= f( x 3~ x i) (y 3 ~yi) > 

(yi~y 2 b 
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APPENDIX D: RESULTS FOR CONDITION NUMBER MINIMUMS 


GPS Tracks and Calculated Track 


for 4 Ranges per Waypoint with DLC Method 



GPS Tracks and Calculated Track 
for 4 Ranges per Waypoint with DLC Method 
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APPENDIX E: CORRECTING EIGENVALUE ERROR 


Actual and Calculated Vnhidn Path lot Simulated Data 



Example of DLC algorithm output after data from matrices with eigenvalues less than 1 
are filtered out. The plots were created using synthetic data with 10 reference nodes. 
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APPENDIX F: STRAIGHT LINE NODE DEPLOYMENT 


Actual and Calculated Vehicle Path for Simulated Data 



Actual and Calculated Vehicle Path for Simulated Data 



When three reference nodes are deployed in approximately a straight line, the FP method 
performs well, whereas the DL method has a much higher average positioning error due 
to inverison of an ill-conditioned matrix. 
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APPENDIX G: PRELIMINARY EXPERIMENTAL PROPOSAL 


Del Monte Lake Initial Experiments APR 2012 

ENS R. E. King, USN 

ESTABLISHING A SEAWEB NETWORK IN DEL MONTE LAKE 

& 

COLLECTING RANGE DATA AND GPS FIXES FOR A SURFACE 

VEHICLE 

OBJECTIVE : To deploy Seaweb nodes in lake and establish a network. To collect range 
data between submerged Seaweb nodes and mobile node, towed by surface vehicle, and 
record GPS fixes along vehicle track. 



[Fig. 1] Sample Arrangement of Submerged Nodes in Del Monte Lake (Image from Google Earth) 

TEST SITE, ENVIRONMENTAL IMPACT & SAFETY : The Del Monte Lake on the 
NPS campus, shown in Appendix (B), is a shallow fresh-water lake, approximately 500 ft 
by 500 ft. The depth is estimated to be between 4 and 7 ft. The lake is fed by drain water 
and has a silt bottom. Historically, this body of water was part of an estuary system until 
it was isolated from the sea by the construction of a railway. The lake was dredged in 
2009. The lake features an island and a fountain that will introduce ambient noise to the 
experiment. The lake is advantageous for Seaweb experiments due to its ready access by 
students, its security as part of the NPS campus, and its lack of marine mammals and 
other species that might otherwise be impacted by acoustic transmissions. Acoustic 
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transmissions are within human hearing range and will not negatively impact the local 
wildlife. The equipment used in this lab is made of non-corrosive materials. Lake tests for 
this experiment are estimated to last one week after which all equipment will be removed 
from the lake. All equipment, including the boat motor, uses clean battery power. All 
provided instructions for motor and battery operation will be followed. Two people will 
be present at all times in the boat. Life vests will be worn at all times. 

READING REFERENCES : 

(a) Endura Pro C2 Trolling Motor Owner’s Manual 

(b) DieHard Automatic Battery Charger Owner’s Manual 

(c) LifeSled Specifications Sheet, www.lifesled.com 

(d) G-Log 760 GPS Recorder User’s Manual, www.transystem.com.tw 

(e) Acoustic Telemetry Modems User’s Manual, P/N 003452, Rev. E 

(f) UDB-9000 Universal Deck Box User’s Manual, P/N M-270-10, Rev. B 

EQUIPMENT : 13 foot Aluminum Boat; Minn Kota ENDURA Pro C2 55 lb Transom- 
Mounted Trolling Freshwater Motor (Model #11910023); DieHard Marine Deep Cycle 
Battery, Group Size 27M (Model #27524); DieHard Automatic Battery Charger (Model 
#71222); LSI towed sled; G-Log 760 GPS Recorders (S/N 8410000916 and 
8410000917); Garmin GPSmap 76CS logger (S/N 10R-022508); Teledyne Benthos 
Transducer (Part #C-270-167-l, Model #AT-440-LF-0, S/N 45444); Teledyne Benthos 
Deck Box (Part #013189 Model # (XDCR) UDB-9000M, S/N 45450); Seaweb repeater 
nodes (Part #013308, Model #ATM-885); 15 lb Greenfield Products Mushroom Weights 
(Part #515-E-UPC); 24 lb Surface Floats. 

Also needed: Oars; Life Preservers; Rope; Battery Cables; Transducer Cables; Assorted 
Tools 



[Fig. 3] The boat used in this experiment is similar to the Sears Gamefisher aluminum boat shown. (Image 

from GovernmentAuctions.org) 

The aluminum boat used for this experiment is borrowed from the Public Works 
Environmental group at NPS. It is approximately 13 ft long and 4 ft wide, intended to 
carry up to 300 lb. 
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Telescoping Handle 



Transom Mount, 
Lever-Lock Bracket 


Flexible Composite 
Shaft, 30" 



Power Prop, 
Weedless Wedge 


[Fig. 4] Minn Kota ENDURA Pro C2 Transom Mounted Trolling Freshwater Motor. (Image from 

www.minnkotamotors.com) 

The trolling motor is a transom-mount design. The depth of the propeller is adjustable. 
The shaft is 30 inches. The shaft is sturdy, but designed to flex in case of impact on rocks 
or other objects in the lake. The propellers are Weedless Wedge 2 type, designed to avoid 
tangling in weeds on the bottom of the lake. The motor is battery-powered and has 55 lb 
of thrust. “Pounds thrust” is the common rating for trolling motors. For a 13 foot boat, 30 
pounds of thrust are required, assuming the boat is carrying two people and little to no 
additional weight. The motor will be used for towing the LifeSled, and the boat will be 
loaded with the weight of the batteries and the modem’s deck box. Operating with 
additional weight, high winds, currents, weedy bottoms, or other adverse conditions 
requires additional thrust. Adverse conditions are not expected on the lake. The 55 lb 
thrust motor should be more than capable of moving the boat and towed equipment. All 
motors will achieve speeds of only 2-4 miles per hour. Trolling motors are designed to be 
quiet so as not to scare fish. The quiet operation will be beneficial for our experiment. 
The motor operates on 12 Volts. It requires one 12-volt deep cycle marine battery. 



Battery Type: 


Group Size: 

27M 

Temperate Zone: 

All climates 

Transportation Type: 

Marine and recreational vehicle 

Power Configuration: 

Voltage: 

12.0 

Power Ratings: 

Amp Hours at 20 Hour Rate: 

105 

Cold Cranking Amps (CCA at 0 deg.F): 

575 

Reserve Capacity (RC): 

180 min. 


[Fig. 5] DieHard Marine Deep Cycle Battery, Group Size 27M. (Image and table from www.sears.com) 
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The batteries used are 12-volt and intended for use with marine vehicles. They each 
weight 57 lb The weight must be taken into consideration for the total weight of the boat. 
Instructions for battery maintenance, operation and safety are included in the DieHard 
Automatic Battery Charger Owner’s Manual (b). 



Amps: 

10/2/50 

Auto Rev. Polarity Protection: 

Yes 

Automatic: 

Yes 

Engine Start Amps: 

50 

Fast Charge Amps: 

10 

Item Type: 

Tools 

Manual Operation: 

Yes 

Overall Color: 

Black 

Primary Type: 

Switch 

Spark Proof Auto Mode: 

Yes 

Trickle Charge Amps: 

2 


[Fig. 6] DieHard Fully Automatic Battery Charger and Engine Starter, for 12-volt Batteries. (Image and 

table from www.sears.com) 

The batteries will be charge on shore over a period of 8-12 hours. The battery charger 
plugs into a standard wall outlet. All instructions and safety measures should be observed 
as outlined in the Owner’s Manual. 



The LifeSled LSI is 62 in x 36 in x 4.5 in and weighs 42 lb. It has over 450 lb of 
deadweight buoyancy. The sled has a 3-point hook up system for towing. It has 10 
internally anchored handles for securing equipment. 


1. Power bottom 

2. Push to log (POI) bottom 

3. Smart mode Led (Blue) 

4. Memory capacity LED (Red) 

5. GPS status LED (Orange) 

6. Car mode LED (Blue) 

7. Battery status LED (Red/Green) 

8. Built-in patch Antenna 

9. USB connector & Reset button (with cover) 


8 


[Fig. 8] LifeSleds 


[Fig. 9] (Image from User’s Manual(c)) 
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The G-Log 760 GPS recorder is a waterproof, all-enclosed device that is 52 x 82 x 19 
mm, and therefore easily mountable to the towed sled. It can log up to 50,000 waypoints 
with a 1 Hz refresh rate and can operate for 25 hours. It uses a rechargeable lithium-ion 
battery and charges via USB connection to a computer. The operating parameters for the 
lake tests, including maximum speed, acceleration, temperature, altitude and humidity, 
are well within the limitations of the logger and the lithium-ion battery. The GPS logger 
records date, time, longitude, latitude, altitude, speed and heading for every fix. It is 
accurate to within 3 meters. 

The G-Log 760 has three modes: normal, smart and car mode. Normal mode records a fix 
at the set interval and remains on until manually powered off. Smart mode adjusts the fix 
frequency based on the current speed and turns off the logger if the logger is stationary 
for 3 minutes to conserve power. Car mode is intended for higher speeds. Data is taken 
every 5 seconds, and only while logger receives external power, usually from the car 
cigarette lighter. Loggers will be used in normal mode to avoid accidentally being 
powered down if stationary during a run, and to ensure a constant fix frequency. 

If logger is in car or smart mode, pushing the reset button (9) will return the logger to 
nonnal mode. NOTE: Pushing and holding the reset button for 5 seconds will return the 
GPS logger to factory default settings and will erase any logged data. Be sure NOT to 
hold reset button if returning the logger to normal mode. 



[Fig. 12] Mushroom Weight [Fig. 13] Surface Float 

(Image from www.westmarine.com) 


The mushroom weights are 15 lb cast iron weights with no-mar plastisol protectant 
coating to avoid corrosion. They are intended for use with sand, clay and firm soil bottom 
types. The lake is known to have a primarily silt bottom. The 15 lb should adequately 
weight down the transducer so that there is no slack in the rope. 

The surface float provides 24 lb buoyancy. The surface floats will serve to keep the rope 
taught and mark the location of the nodes for use when transiting the lake. The length of 
rope used will vary depending on the depth at the site. 
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[Fig. 14] Repeater Nodes 

We are using 7 Teledyne Benthos ATM-885 Acoustic Telemetry Modems (labeled #20-26) 
that each include a built-in omnidirectional transducer. The modem battery pack and 
electronics are contained in a hard coat anodized aluminum housing. The modem has a built-in 
piezoceramic transducer encased in urethane also with a protective anodized aluminum housing. 
The modems can operate at depths up to 2000 m. The transducers operate in the low frequency 
band of 9-14 kHz with a toroidal beam pattern that is near-omnidirectional. The attainable range 
depends on the source level, signal frequency, beam pattern, depth of modem, transducer 
orientation relative to other transducers, and environmental characteristics. Typical horizontal 
ranges are 2-3 km for omnidirectional transducers. The modems can operate independently on 
battery power or can be connected to a host processor and dc power source via a serial 
connection. The modems will be fixed to the anchor weights and surface floats so that they 
are suspended approximately above the lake bottom. 

PROCEDURE : 

Preparing the GPS Recorder for Use - Detailed instructions for G-Log 760 operation, 
software installation and data viewing can be found in the User’s Manual (d). Settings 
and instructions required to complete this lab are detailed below. A summary of LED 
display meanings is included at the end of the experimental procedure document. Refer to 
[Fig. 3] for button and LED locations. 

1. Before leaving the lab, install driver and software for GPS device. With the recorder 
plugged in and powered on, open the TSI LogView program. Click on “Configure 
GPS Module.” Set the fix frequency to 1 fix per second. Do this for each GPS device. 

2. Open the GpsView Program. Baud rate is already set. To find the COM number (for 
Windows), right click My Computer -> Properties -> Hardware -> Device Manager - 
> Ports (COM&LPT). GPS device should be shown with its COM number. 

3. In GpsView Program Status tab, click On. Switch to Setup tab. Fix Update Rate 
should be 1. In AGPS section, click Update. Once complete, you can close the 
program. Do this for each GPS device. 

NOTE: Be sure to turn off the GPS recorder before unplugging the USB cable. 

NOTE: GPS recorders are charged via USB. Ensure both recorders are fully charged 
before leaving the lab. 
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4. Power on the GPS logger by pressing and holding the power button (1) for 3 seconds. 
Press the power to check the battery status when the power is on before each run. 
Orange GPS status LED (5) will stay on while recorded finds satellite signal. This 
may take up to 1.5 minutes. Be sure to power on the device in an open space away 
from buildings with a clear view of the sky. Remain stationary during satellite signal 
acquisition. GPS status LED will blink when signal is found and recorder will begin 
logging. 

5. Push POI button (2) to record Point Of Interest. Memory status LED (4) will blink 3 
times to confirm POI was recorded. If POI is recorded at a kn own location on the lake 
perimeter, GPS data can be compared to the physical location on a Google Earth 
display. Record a POI where the boat is launched and retrieved for every run. Make 
note of physical location for later comparison. 

6. Fix GPS logger to sled. 

NOTE: A third, more accurate GPS unit will be used in addition to the two G-Log 
760 recorders. 


Attaching the LifeSled and Mobile Node - No tools are required for attaching the 
node to the LifeSled other than a crescent wrench to tighten the lock nuts. This can be 
done in the lab, before heading out to the lake. 

The transducer is fixed to the LifeSled as depicted: 



[Fig. 15] Towed sled fixed with acoustic modem 


8. Attach the 25 m cable to the deck box after the sled is sitting in the water. 

Fixing the Motor to the Boat and Attaching Battery Supply - 

9. Fix the motor to the stem of the boat using the simple clamp mechanism included. 

10. When in deep enough water, adjust the motor height to ensure the propeller is under 
the bottom of the boat. 
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NOTE: Be sure to pull the motor up when in shallow water to avoid damage to the 
propeller. 

11. Connect the cable leads from the motor to the exposed connectors on the battery box. 

12. Use the battery box cables to hook up the marine battery inside. The LED indicators 
should light up indicating that the battery is charged. 

Deploying the Seaweb Nodes - Surface vehicle will be operated by two people. Life 
vests should be worn at all times. 

13. Secure one mushroom weight to the end of a rope at least 10 feet long. Mark the rope 
every 6 inches from the bottom of the mushroom weight for at least 10 feet. Use this 
rope for depth measurements. 

14. Load the equipment into the boat, keeping note of the total weight. The boat should 
not carry much over 300 lb total, including the motor, battery, and passengers. 

15. Once the boat is in the desired node location, lower the measurement rope and weight 
into the water until the weight is resting on the bottom of the lake. Record the water 
depth and retrieve the weight. The water depth can also be recorded using a 
handheld fathometer. 

16. Adjust the rope length on the first node based on the local depth. Lower the node into 
the water. Make note of the node number written on the modem. Once the weight 
is settled on the lake bottom, record the node position with all GPS loggers used in 
the experiment. Make visual note of the location of the node for comparison. 

17. Repeat for all seven nodes. 

Operating the Surface Vehicle for Data Collection - Experiment plans will change 

based on the data received form initial tests. The following tests are planned: 


Run # 

Run Description 

Purpose 

Deploy Nodes 

Record GPS log data while 
deploying nodes. 

Ensure GPS recorders are logging 
data as expected. 

1 

Vehicle stationary in center of 
lake for extended period of time. 

Send 4-5 pings and observe 
return data. Repeat for multiple 
positions about lake. 

Initial test to ensure equipment is 
working. Repeated ping and echo 
measurements from same 
location can be compared to 
check that measurements are 
reliable and to get an idea of 
expected error in measurements. 

2 

Navigate vehicle around 
perimeter of lake, navigating 
close to, but not over, nodes. 

Observe what echoes are 
available and reliable at bounds 
of experimental area. Test 
algorithm performance when 
close to nodes. 

3 

60 minute random path with 12 
ping/min 

Collect random path GPS and 
range data for use with Matlab 
algorithms. Collect data that 
spans the entirety of the 
operational area. 


94 









NOTE: Be sure to regularly check battery level for motor, deck box, and GPS units. 


Symbol 

Function / 

Color 

Status 

Description 

i 

Battery Status / 

Red 

Blinking 

Low battery level 

On 

During power charging via USB port 

Battery Status / 
Green 

On / on 

Pow er gauge w hen both POI and Pow er 
buttons are pressed at the same time 

Battery Status / 
Green & Red 

All LED On 

Fully charged 


GPS Status / 

Orange 

On 

Searching for satellite 

Blinking 

GPS fix & start log 

CQ~Q> 

Car mode status / 
Blue 

On 

In car mode, device will go into sleep 
mode w hen no power supply from USB 
and device will wake up and start logging 
when there is power supply from USB. 


Memory capacity / 
Red 

Blinking 

The device is logging position 

E) 

Blink 3 times 

POI button is pressed 

On 

The memory is >80% full 


Keep Hashing 

Erasing data 

© 

Smart Log mode / 
Blue 

On 

Smart log mode on/ vibration sensor on 

Off 

Normal mode on/ vibration sensor off 

2 sec flashing 

Sleep mode 

1/3 sec 
flashing 

Reset button pressed to rest the device 


[Fig. 16] LED Display Summary from G-Log 760 GPS Recorder User’s Manual (d) 
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